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ABSTRACT 

This thesis created an analytical model for active vibration control of the NPS 
space truss using ANSYS. The NPS space truss is a 3.7-meter long truss that simulates a 
space-borne appendage with sensitive equipment at its extremities. With the use of a 
dSPACE data acquisition and processing system, quartz force transducer and 
piezoelectric actuator, he controls using an integral plus double integral control law 
were used to damp out the vibrations caused by a linear proof mass actuator. en 
` reductions on the order of 15-20 dB were obtained with experiment. 

The ANSYS finite element model used SOLIDS elements to model the 
piezoelectric characteristics жа ANSYS Parametric Design Language to provide for an 
iterative approach to an active controls analysis. Comparative data runs were performed 
with the ANSYS model to determine its similarity to experiment. The ETT model 
produced power reductions of 18-22 dB, demonstrating the ability to model the control 
authority with a finite element model. This technique can be used and modified to 
enhance its flexibility to many types of controls and vibration reduction applications. An 


analytical model for active control of the NPS space truss using MATLAB/Simulink was 


also developed as an alternative to the ANSYS model. 
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Е INTRODUCTION 


A. BACKGROUND 


As the 21° Century approaches, the use of space and satellites will continue to 
grow. Even as technology yields smaller devices with reduced power requirements, the 
need for faster and higher-bandwidth communication has led to larger satellites. As 
satellites have become larger and more complex, the launch booster capability has not 
been appreciably increased. This has led engineers to design their satellites to be lighter 
and more compact, and able to deploy large lightweight antennas and reflectors when on 
orbit. The main structure of the International Space Station (ISS), for example, is a 360- 
foot long truss, when measured end to end. To meet the logistics of bringing such a large 
structure to orbit, weight has been a primary consideration. 

These lighter, more flexible structures are prone to low frequency vibration, 
which brings with it new challenges for dynamic control. Accurate modeling of the on 
orbit characteristics of these structures is essential to being able to reduce their vibration 
in the space environment. When on orbit, perturbations may come from a number of 
sources, for example: attitude control maneuvers, crew motion, thermal effects and 
docking and undocking of various orbital transfer vehicles like Soyuz and the Space 
Shuttle. These disturbances may cause large and unwanted vibrations in the structure that 
may disturb sensitive operations. 

In order to 1solate the sensitive equipment from these vibrations, there are three 
possible locations to isolate the disturbance. These options are to isolate at the source of 


the disturbance, isolate equipment at the sink with rack-mounted isolation gear or to 
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remove energy along the disturbance path, usually the spacecraft structure [Ref. 1]. This 
| thesis will discuss vibration isolation along the structural path. 

Inherent in all structures is a degree of natural passive damping. This may be 
enhanced with installed devices, such as visco-elastic dampers. These devices may be 
large, and for larger structures may not be worth the additional lift capacity required to 
bring them to orbit. Active damping is the second option available for vibration 
reduction. 

Implementing an active damping system can be difficult due to the difficulties in 
modeling the dynamic characteristics of the structure. However, once a model is 
developed the control laws and active control devices may then be chosen and evaluated. 

The Naval Postgraduate School (NPS) Space Truss is a small-scale ЈЕ 
structure that has a piezoceramic actuator installed. Previous experiments have shown 
that large reductions in vibration amplitude can be gained with the proper control law 
selection and appropriate gains. Piezoceramic actuators were selected because of their 
high bandwidth and low power consumption. The NPS Space Truss is now a test 
platform for further experiments on active vibration reduction. 

By applying a sophisticated computer-modeling program, a detailed dynamic 
model of the NPS Space Truss can be obtained, and active controls simulated on this 
model. The simulated results can then be verified with an experiment involving the NPS 
Space Truss. The techniques learned in modeling the NPS Space Truss may then be 
applied to larger truss structures, such as the International Space Station. Once a model is 
obtained, a simulation of active controls may be performed, which may yield results that 


are useful for future integration into a structure’s design. 
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B. OBJECTIVES 


The objective of thıs thesis were as follows: 
e Create a finite element model (FEM) of the NPS Space Truss. 
e Develop an active vibration control model and integrate into the FEM. 


e Compare the active FEM to experimental results. 


c SCOPE OF THESIS 


The scope of this thesis includes: 

e Integration of a new dSPACE digital signal processing system [Ref. 2] as the 

heart of the active controls system on the NPS space truss. | 

e Creation and evaluation of an actively controlled а of the NPS space 

truss using ANSYS [Ref. 3]. 

e Comparison of the FEM with the actual truss by experiment. 

The use of ANSYS, MATLAB [Ref. 4] and dSPACE was instrumental to the 
work contained herein. This thesis is not intended to replace the use of tutorials available 
in MATLAB, ANSYS or dSPACE. On-line and printed guides and.tutorials are available 
to give the un-experienced layman an opportunity to learn the system in order to use the 
information presented. A basic outline of how to get the programs and systems operating 
15 included, and has sufficient detail to allow a user with some experience to run the 
desired programs. Where specific program ООО are given, menu commands will 


be used with a “>” to indicate a sub-menu selection, or the commands as typed. 


D. METHODOLOGY 


The research path taken for this thesis followed several parallel paths. The 
research conducted involved the use of active controls systems, FEM creation, fiber-optic 
strain gages, and the specifics of implementation on the programs used. A study of the 
applications of smart structures was conducted in the scope of this research. 

An independent FEM of the NPS space truss was developed using MATLAB to 
verify the selection of the active elements that was installed previously. This FEM was 
also created for integration of an alternate controllable FEM. 

It was decided to create an example active vibration control model in ANSYS 
using generic techniques to test the method and programming algorithm. This example 
model consisted of a single axially loaded strut with a piezoelectric sensor and actuator. 
Once this method of control was verified, it was integrated into a model of the NPS space 
truss. Simulations were conducted in order to evaluate the FEM. 

A new dSPACE digital processing system was installed and tested on the NPS 
space truss to allow for a dedicated hardware/software package. Once the control system 
was installed, it was tested under several experiments to verify its performance, and for 


comparison with the active controlled FEM developed with ANSYS. 


E. THESIS ORGANIZATION 


Chapter II contains background information, theory and a summary of the work 
that preceded the tasks performed by the author. Chapter III 1s a description of the 


installation of the new active controls system, and its implementation, using MATLAB / 


Simulink, and dSPACE. Chapter IV contains the details on the development and 

implementation of an actively controlled FEM of the NPS space truss using ANSYS. 
Chapter V describes the experimental methodology used to venfy the NPS space 

truss with the finite element model created in ANSYS. Finally, Chapter VI will provide 


conclusions and recommendations for further study in the areas researched. 





П. THE NPS SPACE TRUSS 


a TRUSS DESCRIPTION 


1. Background 


The NPS space truss was constructed as a continuation of a series of experiments 
performed by the Naval Research Laboratory (NRL) in Washington, D.C. The NRL 
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space truss was a 3./-meter long, 12-bay squat-“t” shaped structure. It possessed two 
integrated piezoelectric actuators and co-located sensors. The sensors provided input to 
an integral plus double integral force feedback (IDIFF) controller that fed back its signal 
to the actuators. 

Broadband testing was performed from 0 to 100 Hz to determine the operating 
characteristics of the NRL space truss. The engineers at NRL obtained approximately a 
100-fold reduction in amplitude power due to the active control system using the 
controller [Ref. 5]. In cooperation with the Spacecraft Research and Design Center 
(SRDC) at NPS, the NRL provided the main structural components for the assembly of 
an exact replica of the NRL space truss at NPS for further testing. | 

To date, the NPS space truss has had modal testing performed, and a single active 
control piezoelectric strut installed and tested. Current research into the NPS space truss 
‘ involves integration of fiber-optic strain gages and finite element modeling of the truss. 
A catastrophic flood in the adjacent laboratory of the SRDC led to two feet of mud in the 


laboratory containing the NPS space truss. The truss and most of the equipment was 


salvaged, at minimal damage, but most of the connecting cables and documentation was 


destroyed. Part of the work for this thesis ınvolved restoring the NPS space truss to fully 


operational capability. 


2. Elements and Construction 


The NPS space truss is composed of 52 aluminum nodes joined by 161 elements 
in a cubic 12-bay structure. The truss measures 3.67 meters in length and 0.67 meters 
tall. It 1s attached at the center to a fixed base plate. The layout of the unmodified truss is 


shown in F igure 1. The NPS Space Truss. 
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Figure 1. The NPS Space Truss 

The node balls were precision milled from 7075-T6 Aluminum. They are 1.325 
inches across at the threading. There are 18 threaded sites in the three axes and diagonals 
for flexibility in assembly. A standoff was used to provide an interface between the 
individual elements and the node balls. The standoff was fitted with a retaining nut, and 
screwed onto the node ball in the desired location. 

The elements were constructed from 5/16-inch aluminum alloy tubing with 0.035- 
inch wall thickness. They were made in two lengths, a shorter length for the battens and 
longerons (100 total), and a larger length for the diagonal elements (61). A threaded 


sleeve was epoxied to each end of the elements to provide for mounting and precision 


length adjustment. The sleeves were also pinned to the tubes for reinforcement [Ref. 6]. 


An exploded view of the element to node assembly is shown in Figure 2. 


Node Ball (52) 





Figure 2. Node Assembly Details [From Ref. 6] 


During construction of the elements, and during their assembly as a structure, 
great care was taken to ensure that each component was identical. To ensure this, a 
detailed assembly procedure was followed [Ref. 7]. The pieces were assembled such that 
there were no pre-stressed members 1n the truss. This procedure was lengthy, and took 
approximately two man-days to complete. 

Current research at the NRL consists of a second truss of similar "ees 
constructed with carbon-fiber elements. The new truss was constructed without the 
previous rigorous torque specifications. This would better simulate the on-orbit assembly 
by an astronaut performing extra-vehicular activity (EVA). For an astronaut performing 
an EVA, some, but not all, of the assembly is rigidly specified. 

After assembly, the four lower-central nodes were rigidly attached to a base plate, 
for mounting on 4 Newport Vibration Control System Table. The 1000-Ib. isolation table 
uses compressed Nitrogen to charge pistons in its legs, providing a cushion under any 


mounted equipment and providing for high frequency vibration attenuation (greater than 


О) 


99% above 12 Hz) [Ref. 8]. Detailed data on the NPS space truss are contained in 


Appendix A. 


B. ACTIVE CONTROL ELEMENT 


1. Overview 


The active control element for the NPS space truss consists of a piezoelectric 
actuator and co-located force transducer. These components were mounted to steel rods 
in a manner to be compatible with the truss elements. Together with the computer data 


processing system, these components comprise the active control system. 


2. Piezoelectric Theory 


Piezoelectric (also called piezoceramic) materials are materials that will elongate 
when an electric field is applied in a pre-determined directión. Conversely, when a 
deformation is applied to a piezoelectric material with either an external force or strain, 
an internal electric potential will be developed that can be measured. In this manner the 
piezoelectric material can be used both a structural sensor and as an actuator. 

Some piezoelectric materials, such as quartz, occur ПЕШІ) others must be 
synthesized [Ref. 9]. A piezoelectric material is created when a suitable material, such as 
Lead Zirconate Titanate (PZT) is heated above its Curie point under an externally applied 
electric field. The Curie point is the temperature where atomic magnetic alignment 
becomes unfixed. An applied electric field at this point will permanently realign the 


spins and magnetic dipoles [Ref. 10]. Cooling while maintaining the field will lock the 
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magnetic dipoles in place. This creates the poling axis, P, where elongation will occur 


under an externally applied electric field, E. 

Piezoelectric materials are anisotropic in nature. The directions that are 
orthogonal to the poling axis are related mechanically by Poisson’s ratio, in a similar 
manner to purely mechanical deformations. Applied deformations in the directions 
orthogonal to the poling axis will cause an electric field to be created in the poling axis 
due to the mechanical coupling. 

While not truly linear, and containing some hysteresis, the piezoelectric effect can 
be modeled as such, with little error. The piezoelectric electroelastic relations can be 


described with the following equation [Ref. 8], shown in stiffness method form: 


0) id ШЕ di 


The stress and strain vectors represent the three axial and three shear stresses and strains. 


Equation (2.1) can also be represented in an applied force method, given by [Ref. 11]: 


foliar 219 aa 


These equations, when the [d] and [e] matrices are zero, reduce to Hooke’s Law and the 

dielectric equation [Ref. 9, p. 4]. The relationship between the [d] and [e] matrix is given 
by the following: 

[e] = [c] [d] (2.3) 

The dimensions of the [c] matrix are six by six. It contains all of the linear elastic 


and flexural terms, and is symmetric. The [e] and [d] matrices are six by three in 
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dimension and contain the piezoelectric coupling terms. Finally, the [e] matrix is a 
diagonal matrix with the dielectric permitivities. The individual constants in the matrices 


are usually indexed by their coordinate axes (1, 2 and 3); the 3-axis is the poling axis. 


3. Application and Operating Characteristics 


Piezoelectric sensors and actuators are attractive for use for a number of reasons. 
As a strain gage, a piezoelectric material can be as high as one millton times more 
sensitive than a traditional metal-foil strain gage. They also possess low noise, and being 
a ceramic material, low to moderate temperature sensitivity. Their bandwidth and 
response are compatible with their use as actuators. The power required by a 
piezoelectric actuator is low, but due to the high voltages required, requires special 
attention when placed in an orbiting spacecraft. 

Basic schematics for a piezoelectric sensor and actuator are presented in Figures 3 
and 4. The sensors and actuators can be mounted such that the sensing direction and 
desired applied force can be in the direction of the poling axis, or orthogonal to it. In the 
following figures, the poling axis and the applied or sensed electric field are shown, as 
are the displacements in the piezoelectric crystal. 


-У 





Figure 3. Piezoelectric Actuator Configuration 
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Figure 4. Piezoelectric Sensor Configuration 
Piezoelectric materials are limited in that the electric field that can be applied 
must be less than 2kV/mm [Ref. 12]. To functionally eliminate this problem, several 
sans of piezoelectric material are stacked in the direction of the poling axis, with 
interleaving foil electrodes to provide a greater cumulative strain and thereby greater 
force to the structure. The polarity of the piezoelectric wafers is inverted at each step to | 
simplify the electrode placement and maintain uniform global effect. This configuration 


for a stacked linear actuator 1s shown in Figure 5. 





<= Defines polarization direction in wafer 


Figure 5. Stacked Piezoelectric Actuator [From Ref. 13] 
The active piezoelectric element is mounted in the structure such that the poling 


axis 1s directly in line with the element that it is replacing. This mounting will allow the 
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maximum piezoelectric effect to be felt by the structure, and minimize damage to the 
brittle ceramic wafers within. 

The piezoelectric active control element that was selected was the Piezoelectric 
Translation Model P-848-30, built by Polytek Physik Instrumente of Hamburg, Germany. 
The P-843.30 has a maximum operating voltage of 100V. It is cylindrical in shape, with 
a 14-mm diameter su 73-mm length. The P-843.30 can be seen in Figure 6. Other 


characteristics are contained within Appendix A. 





The P-843.30 has two connectors, one for the piezoceramic wafers, the other for 
an integrated metal foil strain gage that can be also used for controls applications. The 
strain gage was not used for these experiments. A Fitting was added to make the 
connectors compatible with a BNC connection. | 

The P-843.30 has a rated open loop travel characteristic of 45 um/100V +/- 20- 
percent. An experiment was performed to i this amount for the actuator prior to 
installation [Ref. 14, pp. 25-28]. To perform this test, the actuator was mounted to a 
rıght-angle stand and a voltage applıed. The dısplacement was measured wıth a Kaman 


Eddy Current Sensor. The data recorded was within the manufacturer’s stated tolerance 
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and displays the hysteresis discussed earlier. The results from this measurement are 
presented in Figure 7. Piezo Model P-843.30 Expansion Characteristics [From Ref. 14]. 
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Figure 7. Piezo Model P-843.30 Expansion Characteristics [From Ref. 14] 
As can be seen from the previous Figure, there are two piezoelectric actuators. 
One of the actuators was damaged irreparably during installation; the serial number on 


the installed actuator can not be read. It was assumed that the average value would 


suffice for further calculations, at 50 um/100V. 


4. Locating The Active Element - 

A finite element model of the NPS space truss was constructed in MATLAB to 
verify the active strut location. It was created using the three-dimensional frame elements 
found in Kwon and Bang [Ref. 15, p. 264]. These elements are subjected to axial, 


bending and torsional loads. For a complex shape, like an I-beam, the elements can be 
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complex. In the NPS space truss, the elements are all tubular and axially symmetric, and 
some simplifications were used. 
The equation of motion for a simple n-degree of freedom undamped system is 
given by the following equation, shown in matrix tensor notation: 
[Mi +[K](x) =(F) (2.4) 


The stıffness matrıx, [K], for a three-dimensional tubular frame element 1s given 





by 
ESIE В An | (2.5) 
К», К», 
where 
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The constants in the stiffness sub-matrices reflect the different effects of 
compression, torsion and bending. The ‘a’ terms reflect the axial stiffnesses both 
compressive and torsional. The ‘b’ and ‘c’ terms are used for the bending properties of 
the frame elements. The three ‘b’ and ‘c’ terms reflect the superposed sum of deflections 
from two displaced nodes with zero slope, one displaced node with zero slope and the 
other with slope but no deflection, and zero deflection but with slope on both ends. 


These terms are all described below: 


y AR ios NN boe T. (2.7) 





For each element, there are twelve degrees of freedom: six that relate to the three- 
axis displacements at the two nodes, and six that describe the rotation of the nodes in the 
three axes. The shape functions for the axial displacement and rotation are linear, while 
the off-axial displacements and rotations are cubic, Hermitian polynomials, which are 
independent from each other. 

As each elemental stiffness matrix 1s formed, it is 1n a local coordinate frame, with 
the x-axis ТШ with the long axis of the element. The y- and z-directions are 
orthogonal to this in a user-defined direction. The other axes are important when the 
beam is not symmetrical about the axial direction. To rotate the coordinate frame from 
the local to the global requires the use of a direction cosine matnx (DCM). The DCM 
consists of the dot products from one coordinate frame to another, in this case the (u, v, 
and w) local frame and the (1, j, and k) global frame. This rotation can be seen in matrix 


tensor notation below 
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[к =[Dcm]"[k*] [Dcm]* (2.8) 


For each node’s displacements and rotations, there must be a DCM. The overall 


DCM then becomes 
[аст] 
аст}? 
[рсмү = Рет (2.9) 
gl 
[dem | 
| [аст 
where each individual direction cosine matrıx (dem) is given by 
i 7-08 1.58 1-8 | 
[dem] =| jü 7-08 јен |. (2.10) 
5.08 К-ў* Еу | 


The local coordinate system is defined as having the x-axis aligned with the long ` 
axis of the element, with the y- and z- directions orthogonal to it. When writing an 
algorithm to automate this for a computer program, the 0, v and w unit vectors were 
defined in terms of the global frame. This allowed the cosines (dot products 1n the [dem] 


above) to be easily calculated. In the global frame, the axial direction was defined as 
ыма ые een (2.11) 


Since the beams were symmetric, the direction of the y- and z-local directions was 
arbitrary. This allowed the use of an off-axis vector, which was used іп forming the other 
orthogonal axes, to be chosen by a simple rotation of the x-local direction vector by 90- 


degrees along the z-global axis. This was useful except when the x-local axis was already. 


aligned with the z-global direction, in which case the x-local axis was rotated in the y- 


global direction. The off-axis vector was defined by 


0 1 0 (Жош. 
ОА=|-1 0 0 |28, ог|0 1 0124. (2.12) 
001 100 


The v-direction was derived by obtaining the cross product of the u-vector with 
the off-axis vector. Similarly, the w-vector was obtained by completing the nght-hand 


rule involving the u- and v-direction vectors. These relations are 


ae > 
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With the element DCM obtained, the local stiffness matrix may be rotated into the 
global frame for integration into the global stiffness matrix by using the tensor equation 
(Equation (2.8)). 

The next requirement for structural dynamic analysis 1s the mass matrix, [M]. 
There are many approximations available to a modeler for a mass matrix; some of them 
are lumped, consistent and compound mass matrices. For the NPS Space Truss FEM, all 
three were EA 

The lumped mass matrix consists of taking the total mass of the element and 
assigning half of the mass to each node in each translation degree of freedom. The inertia 
terms are neglected for a lumped mass matrix analysis. The em form of the 


lumped mass matrix 1s given by: 
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1 0 0 0 

l ALIO 0 00 
[m] = (2.14) 

2 Оо 0 TE 

0 0 0 0 


In a multiple degree of freedom system, the matrix is expanded so that each 
translational degree of freedom is given one half the mass of the member. The total mass 
in all the purely translational terms is equal to the total mass of the element, multiplied by 
the number of translational degrees of freedom per node. This is true for all mass 
matrices, regardless of construction method. 

Integrating the shape functions for each degree of freedom, multiplied by the 
mass, over the length of the element forms the consistent mass matrix. The matrix tensor 


form of the consistent mass matrix 1s given by 
[Me] = [ epu] [5 |а, (QUIE 


where [H] 1s the matrix containing the shape functions for the element, applied for the six 
degrees of freedom. The shape function matrix 1s of size one by twelve for a two-node, 
six-degree of freedom element. For the linear axial shape functions and Hermitian off- 


axis shape functions, the [H] matrix 1s in the form 


ІНіІ-ін н н, н Нн HH, Нн AS 


ІС 1С IL 2С 2C 21. зс эс 21, 


. where the linear and cubic shape functions are those found in Ref. 15. The independence 
of the axial displacement, axial rotation and'off-axis displacement/rotations 1s apparent 


from the reuse of shape functions in equation (2.16). 
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It was found through trial and error that including the small rotational inertia 
cross-terms 1n the consistent mass matrix led to floating point errors during the MATLAB 
model run during eigenvalue solving. Therefore, it was decided to neglect these terms 
from further analyses. The form of the consistent mass matrix used in the NPS space 


truss model 1s given by: 


00 0 9990800735. 0 0 000) 
078 0000027 0 000 
Qm cc i 
ОООО ИИ 
000000000000 
[m] O II о ONUS OO] (2.17) 
Оооо ооо ОО 
пов” 
В ооо а” 
ТУО О ООО ОГО ОО 
Мм о о с о о о о 
0770970 0 90070 077070 07% 


The differences found in the axial (diagonal) and off axial (off-diagonal) terms are 
due to the fact that the axial shape functions were linear, and the other directional shape 
functions were cubic. For matrix inversion, the zero rows and columns in the mass 
matrix were removed, and subsequently replaced to restore the degrees of freedom. 

Some re analysıs programs, such as MSC/NASTRAN contain an option 
for a different kind of mass matrix, called a compound mass matrix [Ref. 16]. Taking the 
average of the lumped and consistent translational terms forms the compound mass 


matrix. The form of the compound matrix for linear elements is given by 


туташ 

[м“] =^&|° 0. e (2.18) 
ML D SE 
0 0 0 0 
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For cubic shape functions, the form of the compound matrix is different. As an 
approximation, and to ease computational burden, the compound mass matrix shown 
above was finally used for the NPS space truss analytical model. Finally, the additional 
mass from the node balls and fittings was added as concentrated mass, located at the 
nodes. 

The damping effects in the NPS truss were neglected for the MATLAB FEM, as 
they were not expected to contnbute heavily either to the natural frequency or the best 
location for the control strut. In previous experiments, the damping ratios were found to 
be less than ten-percent, which would lead to less than one-percent error in the natural 


frequency. This relationship is shown by the following equation: 


ЕЕ 0.99 = 0. (218) 

For the finite element modeling, the MATLAB m-file “ex895.m” [Ref. 15, p. 276] 

was extensively modified to form the backbone of the model creation as 

"Final root nps.m." This program is included as Appendix B. In order to automate the 

calculation of the elemental mass and stiffness matrices, ‘feframe2.m’ [Ref. 15, p. 278] 

was modified to provide a three-dimensional element, convert it to the global coordinate 

frame and provide the mass matrix for each element. This modification is found as 
‘feframe3.m,’ included as Appendix C. 

When the root program was run, the eigenvalues of the system were determined. 

The first ten modes are shown in Table 1, NPS Space Truss Modal Frequencies. For 


comparison, the natural frequencies for the actual truss, as determined from a modal 


~ 
~ 


testing experiment [Ref. 14] are presented alongside the frequencies obtained in the 


analytıcal model. 
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Table 1. NPS Space Truss Modal Frequencies 







As can be seen from the table above, there 1s a small difference between the ual 
and calculated natural frequencies. There are several possibilities as to where the 
difference 1s generated. The modeling of the truss as a series of single elements between 
the nodes 1s approximate, and does not include any stiffness for the fittings that join the 
elements to the Aedes The lumping of mass at the nodes also affected the results. In 
later, more detailed models, additional stiffness elements could be added for the fitting 
between the element and the nodes with more model nodes at these fitting locations. 

In any true structure, there are infinite degrees of freedom, but this is analytically 
impossible to model, even for the largest computer. For this model, the results obtained 
are very close. This is likely due to the high number of degrees of freedom of the system 
(312), which enables a better representation of the early modes. 

The first four mode shapes for the NPS Space Truss were obtained by taking the 
eigenvector displacements and adding them to their respective nodal coordinates. This 
will then provide a figure that HOw the position of each node in the given mode shape. 
Each mode shape was created using the code contained in Appendix D, “NPS_modes.m.” 
These figures are included in Appendix E. Modes one and two represent a side-to-side 


and “see-saw” motion of the truss respectively. Modes three and four represent an “arms 


p 


waving” description, sideways and up and down, respectively. 

In order to determine the best location for an active control strut, the following 
approach was used. The structure was analyzed for the first five modes to determine 
which element had the maximum strain energy for that mode. To target the first mode, 
the element with the maximum strain energy should be replaced with an active member. 

When a structure 1s vibrating in any of its various modes, each element has some 
strain energy, caused by displacement along the element length. This energy is at a 
maximum when the structure is at its largest amplitude. The elemental strain ШЕ 15 


given by the following equation: 


— 


Ce. -—E e’= leds: > (2.20) 
2 DL 
where 
6“ -(х-х) +(у,–у) +(2,—4) (2.219) 


The change in element length was determined by using the nodal coordinates, 
when displaced, from the mode shapes determined above. Since eigenvectors are 
amplitudes for ee functions, the displacements can be negative as well as positive. 
This means that the geometrically symmetric elements are also as effective as the truss 
vibrates in the other direction through a cycle. The numerical results for the relative 
strain energy determination were generated using code contained in Appendix F. 

Due to the diagonal elements in the truss, there 1s an odd symmetry to the 
locations of the maximum strain energy. The highest elements occur near the center of 
the truss in all modes except the third (torsional) mode. With this information the 
location for the active control strut may be selected. Element 101 is the location with the 
maximum strain energy for mode one. A second piezoelectric control strut can be placed 
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on the opposing side of element 107. This is possible due to the symmetry described 


above and the desire to not have both actuators acting through the same node ball. 


5. Force Transducer 


A large issue for any active structure’s design is the location of the sensors and 
the actuators [Ref. 15, pp. 460-461]. The vibrational motion of a structure involves 
phasing difference that can vary for all frequencies that a structure will experience. For 
this reason, a co-located sensor is attractive, as it will read the same phase as the applied 
force. In a noiseless system, a co-located sensor can mathematically provide 
unconditional stability, given proper control law. It also provides for analytical 
simplifications due to the application at the same point as the sensed output. 

The force transducer that was selected was the PCB Piezotronics Model 208B02 
General Purpose ICP Force Sensor. This sensor is a linear, piezoelectric device, with 
quartz as the piezoelectric material. It is threaded on either end so that it can be placed in 
the structure to experience an axial load. It is capable of detecting a force of up to 1000- 
lb compressive under static conditions and 100-Ib under dynamic conditions [Ref. 17]. 
The device is 15.75-mm long and approximately 17-mm in diameter. The model 


208B02’s operating characteristics are included in Appendix A. 


6. Installation of the Active Element [Ref. 14] 


Once the active strut components were gathered and tested, they were assembled 
into a single active truss element. Machined steel rods were created to provide an 


interface between the active element and the truss. Given the susceptibility to moment 
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and torsional loads, a PI flexible tip was inserted between the sensor and actuator. The 


assembled active element is shown in Figure 8. 
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Figure 8. Assembled Act 
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To utilize the full range of the piezoelectric device, a preload was required. The 
options were to mechanically preload the device, at the cost of loading the structure, or to 
determine a means for an electrical preload. Mechanical preloading.was used and was 
achieved by inserting shims into the element to node ball interface until the desired 
preload was achieved. Full details on the installation of the active control element can be 
found ın Reference 14. 

During installation of a second active control strut, excess torsıonal load was 
applied, resulting ın failure of the second strut. In the future, care must be take to ensure 
that torsıonal loads are balanced during assembly. The actuator is machined to allow the 


use of an open-ended wrench for this purpose. 
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С. LINEAR PROOF MASS ACTUATOR (LPACT) [Ref. 14] 


In order to provide a disturbance force, a Linear Proof Mass Actuator (LPACT) 
was installed in the NPS space truss. The LPACT ıs a model CML-030-020-1 
manufactured by Planning Systems, Incorporated, from Melbourne, FL. The LPACT is 
powered by a separate amplifier and controller assembly that has embedded feedback 
electronics. These feedback loops, while not used for these experiments, could also be 
used for active damping control. The LPACT also has mounted accelerometers for use in 
monitoring or driving the feedback loops. Details on the LPACT are contained within 


Appendix A. The LPACT is shown in Figure 9. 
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The LPACT was mounted in a similar manner as the active strut. The device was 
balanced to locate its center of mass, to ensure that symmetry would be maintained in the 
truss and to prevent a torsional static load down the length of the truss. Aluminum 
interface struts were constructed and attached to the LPACT for mounting. It was placed 


in the truss in the outermost diagonal, to simulate a payload that is vibrating at critical 


D 


frequencies. After mounting, the Gravity Offload Spring was used to retum the LPACT 
assembly to the center position, to offset the weight of the mass at a 45-degree angle. A 


picture of the installed LPACT is provided as Figure 10. 
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From Ref. 14] 

The LPACT, when the internal feedback loops are disabled, is powered by an 
external signal, from a signal generator, or other source. As the LPACT is a real device, 
it has its own frequency characteristics. A description of the output force and phase are 
shown in Figure 11. It is of note that the LPACT has its natural frequency at about 8.5- 
Hz, and a tapering off of the amplitude from that point. For the experiments to follow, 
the range of 10 to 20-Hz is of particular interest. The mass of the LPACT does contribute 


to a shift in natural frequencies of the truss. This is in addition to the inclusion of a new 


natural frequency in the 8-10 Hz range that is the product of the LPACT spring. 
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Figure 11. LPACT Transfer Function (feedback loops off) [From Ref. 18] 
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D. LASER DIODE ASSEMBLY [Ref 14] 


To allow a qualitative measurement of the truss' vibration when under the 


influence of the LPACT, a laser diode was mounted to the end of node 52, such that it 


would project a beam onto the wall. The mounting consisted of a thin aluminum rod, 


with a mass at the end, housing the laser. The laser diode selected was a 1-mW, 635-nm 
Model PLC6351FW from Lasermate Corporation of Walnut, CA. It 1s powered by a 


Hewlett Packard E3615A DC Power Supply, operating at 2-4 Volts nominal. A picture 


of the laser diode assembly is included as Figure 12. Full details on the laser diode can 


be found in Reference 14. 
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E. PREVIOUS EXPERIMENTS 


Two previous theses have been completed on the NPS space truss: One by CAPT 
Brent K. Andberg, USMC, and the other jointly by LT Scott Johnson and LT John 


V]attas, USN. 


1. Andberg [Ref. 19] 


In his thesis, CAPT Andberg developed an FEM of the NPS space truss using an 
NRL code entitled NRLFEMI. He then performed modal testing on the truss to confirm 
the model. The experimental data was lacking in that it failed to observe the first mode. 
Finally, CAPT Andberg performed a technology demonstration of the use of Fiber-optic 
Bragg Gratings (FBG), used in this example to detect the motion of a simple cantilever 
beam. In the future, FBGs will be installed on the NPS space truss for shape 


determination. 
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2. Vlattas and Johnson [Ref. 14] 


This Master’s Thesis focussed on two areas, the re-performance of modal testmg 
with the HP 35665A two-channel spectrum analyzer and integration and testing of an 
active control device. They recommended that modal testing be again done, to overcome 
some of the limitations of the dSPACE system that was installed at the time. As will be 
discussed in the next chapter, a new dSPACE system was installed on the NPS space 
truss, and would be available for this purpose, if desired. 

The active controls integration produced good results, with a maximum ee 
reduction of 14.817 dB at 16.85 Hz. They reported using a disturbance amplitude of 100 
mV for the LPACT source. It was determine early on through the current course of 
research that this did not even provide а or signal sufficient to overcome system 
noise. Therefore, the 14.817-dB reduction 5 amplitude is held in question. One of the 
experiments performed later will be to validate these results with the new dSPACE 


system. 
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III. CONTROLS SYSTEM INTEGRATION 


A. OVERVIEW 


The active vibration controls system was created with the use of a dSPACE DS 
Оре Controller running the control law ımplemented by Vlattas and Johnson [Ref. 
14]. The DS 1103 has a Motorola PowerPC 604e microprocessor as ıts central 
processing unit (CPU), and resides in a triple-wide ISA slot in a host PC [Ref. 2, p. 12]. 
The dSPACE system described in References 14 and 19 was originally shared er 
several systems; the DS1103 system was ordered and installed for the sole use of the NPS 
Space truss. | 

The DS 1103 board has connectors for an external input / output (I/O) box that 
contains BNC fittings and standard computer cable connections for analog to digital 
(ADC), digital to analog (DAC) and support for other cabling formats (e.g. RS-232). 
Control of the dSPACE CPU and access to its memory (128MB) is done with the use of 
ControlDesk, the dSPACE main program, or by using MATLAB / SIMULINK programs 


or C-code. These programming components are discussed in grater detail below. 


B. INTEGRAL PLUS DOUBLE INTEGRAL CONTROLLER 


The selection of the integral plus double integral force feedback (IDIFF) 
controller was due to its inherent stability. A simplified block diagram of an integral 
controller acting on a unity system is shown in Figure 13. A plot of the magnitude versus 


frequency for this system 1s presented in Figure 14. 
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Transfer Fcn 


Figure 13. Integral Feedback Controller (Simplified) [After Ref. 5, p. 8] 
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Figure 14. Effect of Integral Gain on Closed Loop Transfer Function [After Ref. 5, p. 8] 


In theory, the controller gain would be increasingly able to control the system 
output. However, the authority of the actuator (how much force ıt can ımpart to the 
structure) and the noise in the system lead to a point where the maximum gain is achieved 
[Ref. 5, p. 8]. Further, the addition of digital filters to the design to attenuate the DC bias 
will add additional factor that also invalidate the unbounded controllability inherent in the 
. integral controller. The use of the second integrator can be used to finely adjust the phase 
response of the controller, and allow a slightly better response than the single integral 


alone. This result will be discussed in more detail later in the thesis. 


34 


С; SIMULINK /RTI 


1. Simulink Controller 


The first step in creating the controller was to assemble a Simulink model. The 
Simulink model contains the processing connections necessary, implemented with block 
diagrams, to create the control law. The original controller, found in Reference 14, was 
recovered and updated for the current version of MATLAB and Simulink. 

The original controller contained two band-pass filters designed to isolate a 
desired frequency (in their case, 16.75 Hz). It was decided to test the controller over a 
larger range of frequencies, therefore the band-pass filters were at first removed. It was 
found through several test runs that the input and output, which will be fully described in 
the next subsection, contained a DC bias that needed to be removed, lest the integrators 
give an unbounded output. Therefore, it was decided to implement a high-pass filter in 
the signal path. This need is consistent with the results observed by the researchers at 
NRL, who used a second-order Butterworth filter to remove the DC bias ЕЛ о |. 
Higher order filters were tested to determine their utility. It was discovered that the 
higher-order did not lead to an appreciable increase in filtering, but did decrease the 
instantaneous response time. A third-order filter was used for the NPS space truss. 

The Simulink model also contains a provision to prevent supplying an over 
voltage to the piezoelectric strut. A saturation limiter with DC bias is included as part of 
the control signal, before it reaches the output block. 

The dSPACE system relies upon the Simulink model for the selection of its 


sampling frequency. It was determined that frequencies above 250 Hz (as reported in 


ВЭ 


Chapter V, section B.5) were not consequential. Therefore, a sampling frequency of 
500Hz was utilized, yielding a sample time of 0.002 seconds. When the Simulink model 
is complete, the command “Tools > RTW Build” will create the real-time program and 


object file that can be executed by the CPU of the SPACE system. 


2. RTILIB 

The interface from Simulink to dSPACE is found in the Real-Time Integration 
Simulink Library (RTILIB). It 1s invoked by typing “rtilib” Ота the MATLAB 
command window after Simulink 1s running [Ref. 20, p. 17]. The RTILTB contains all 
the additional blocks required to provide the interface between the Simulink model, 
running on the CPU, and the external I/O box. It also provided blocks for hardware 
interrupts for more complex programming needs. The RTILIB interface is shown in 


Figure TS" 
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To use the blocks ın a model, they are selected and then “dragged and dropped” 
onto the Simulink model for integration. The I/O box connections are found under the 
“MASTER PPC” icon. There are twenty available inputs and eight outputs available to 


the userin dSPACE. They are accessible by three blocks, shown in Figure 16. 


mne mnt e note tnr I i e e t MH t e I MM 





51103 РРС Controller Bord 
Master РЕСНИК 


Pe м ~ ADAN a еі Коми 
ye Ep om уч ттс 


АРС 1 MUXADC pb DAC 


 DS1103ADC C17 DS1103MUX ADC CON1' DS1103DAC C1. 


enims М... р.а o D etg cuts REP па лаке imde aa D 


Figure 16. MASTER PPC VO Connections [After Ref. 20] 

The ADC channels available are the ADC and MUX (Multiplexed) ADC, shown 
above. The MUX ADC is a series of four channels, with four multiplexed inputs each 
(ADCH 1-16). They can be used in Simulink as individual separated channels, or as a 
vectored input. The scalar ADC block is a single input, corresponding to the last four 
input channels (ACDH 17-20). In similar manner, the eight output channels are accessed 
with the DAC blocks (DACH 1-8) [Ref. 20, p. 24]. The I/O blocks, when connected to 
the I/O interface box, exhibit a 10:1 gain in value. For example, a 0.5 scalar value in the 
running model will produce a 5.0-Volt output at the DAC. 

Additional gain blocks were added to the model to allow for turning on and off 
the controller without modifying the gain settings. Also, vectored input blocks were 
added to allow the capture of the four accelerometers used in the experiments. The final 


controller block diagram is shown in Figure 17. NPS Space Truss Controller. 
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Figure 17. NPS Space Truss Controller 


D. dSPACE CONTROL DESK 


In previous versions of dSPACE, data input and output was managed with the use 
of two different software packages: Trace and Cockpit. m. allowed for real-time 
detection of the signal and E variables with software capture, Cockpit А quası 
real-time adjustment of the parameters running on the CPU. In the latest versions of 
dSPACE, these programs have been melded into a single entity, called ControlDesk, and 
expanded in scope to include some of the Windows-type features that the majority of 
computer users have become familiar. In addition to allowing the access and variation of 
data in the CPU, ControlDesk also allows for grouping relevant files under an 
“experiment” and running of macros to automate the data taking process using the 
dSPACE macro language, Python. [Ref. 21, pp. 13-14] 

ControlDesk’s most useful improvement is the addition of Windows “drag and 
drop” (D&D) capability. It allows for the selection of variables for real-time display in 
on-screen instruments, and for entire programs to be downloaded and executed by the 
CPU [Ref. 21, p. 31]. To establish a working experiment, first a program must be loaded 
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into the CPU (or the CPU’s memory), then an instrument panel may be built that will 
monitor and control the parameters and signals desired. The ControlDesk interface 
contains several sub-windows that allow a user to run his/her experiments. These 
windows are the main window, Navigator, and Tool Window. The Navigator allows file 
and program manipulation; the Tool window allows the selection and control of the 
various variables in a model. The main window Alene for the controlling of the 
simulation (i.e. running the CPU) [Ref. 21, pp.43-50]. 

Instrument panels are built using a pull-down menu interface and using D&D to 
place the desired instrument on the main window. There are many different types of 
instruments available, to simplify the visual arrangement of the experiment to the user. 
Some of the instruments available are sliders, pointers, oscilloscopes, knobs and buttons. 
Selecting the proper variable is important. The variables can be of many different types, 
and are grouped according to their block location in the onginal Simulink model. The 


relevant types to this thesis are shown in Table 2. 


| B: | Block Outputs | 
|. P: | Block Parameters | 


Table 2. dSPACE Variable Types [After Ref. 22] 








The controller layout used for the IDIFF controller on the NPS Space Truss is 
shown in Figure 18. 

The four large plots on the top of the Layout are oscilloscope plots of the four 
accelerometers, with the three axes measured superimposed on each plot. The plot in the 
lower right-had portion of the layout is used to track the incoming signal from the force 


sensor and the signal to be output to the actuator. The lower left-hand portion of the 
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layout contains controls for the various gains and a read or green light to indicate from a 


glance the controller status (controlling or uncontrolled). 
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Figure 18. NPS Space Truss Controller Layout 
To start and run the controller, the following steps must be followed (all on the 


NPS space truss computer in the Smart Structures Laboratory): 


е Initiate ControlDesk using the 2 | icon either on the desktop of the toolbar, or 
by using the menu Start>Programs>dSPACE>ControlDesk 
e Open the Experiment using the command File > Open Experiment. Select 


C:\Space_Truss\truss99\Space Truss.cdx. 
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me 


e Load the controller to the CPU by using the = icon on the toolbar. Select 
the program C:\Space_Truss\truss99\ST_Controller_3.ppc (the number refers 
to the third iteration of the controller). 


е Ensure the Edit Mode is selected from the toolbar (verify the Edit mode button 
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depressed) — . The three buttons are Edit, Test and Animate mode 
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respectively. 
e Start the CPU (depress the green triangle on the toolbar). 
e Select the Animate mode to enable the display. 
e The controller may be turned on and off with the push buttons in the lower left 


corner of the layout. 


E. AMPLIFIERS AND CONDITIONERS 


The final segment of the NPS space truss controller consists of the amplifiers and 
signal conditioners that are external to the PC and the dSPACE system. These 
components are described below. The signal that is produced by the force transducer is 
fed into PCB Piezotronics Model 484B signal conditioner, nen has a unity gaın and a 
selectable bias (either 6.0 or 11.0 Volts, DC or AC). The four accelerometers used have 
three-axis capability. Each axis from the four accelerometers is fed into a 12-channel 
Kistler Piezotron Coupler, Model 5124A. Finally, the control signal is sent to a Trek 
50/750 Voltage Amplifier. Details of the equipment are included in Appendix A. NPS 


Space Truss Characteristics and components. 
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The Trek amplifier was due for calibration in September 1998. Due to the 
timeliness involved, it was decided to not perform a full calibration on the Trek amplifier 
and instead perform a mini-calibration on it to verify its linearity and set points. The 
Trek was calibrated for a 20-time voltage gain. The gain verification was done by 
supplying a voltage from 1.0 to 6.0 V to the amplifier and measuring the output voltage. 


The details of the mini-calibration are presented in Table 3. 


0.998 V 20.26 V 
2.006 V 41.2 V 


Table 3. TREK Voltage Amplifier Mini-Calibration 





The deviation from linearity in the data obtained was deemed to be minimal; therefore the 
mini-calibration was determined to be sufficient to allow use of the Trek Voltage 


amplifier for the experiments. 


IV. CONTROLLED FINITE ELEMENT MODEL OF THE NPS SPACE 
TRUSS 


A. OVERVIEW OF ANSYS 


ANSYS ıs a powerful finite element utility that has many capabilities. The 
program can involve many different physical properties, from over 100 types of elements. 
In this manner, the user can model systems that involve, for example, mechanical, 
electromagnetic, thermal, piezoelectric and electric characteristics. These options 
available are based on the licensing obtained from ANSYS. ANSYS can couple these 
effects into one single model, and provide insight into it inner workings. For the work in 
this thesis ANSYS version 5.5.2, Multiphysics was used. 

ANSYS may be run from a graphic user interface (GUI) or purely by text 
commands. In the GUI, there are several menus and toolbars that are present to allow a 
user to perform a number of tasks. These menus and toolbars are as follows: 

e Utility Menu: Allows for the control and manipulation of variables, files and 

graphic displays. 

e Input Window: An interface for typing in text commends directly while in the 

GUI. 

e Toolbar: Contains standard and user-programmed macro buttons. 

e Graphics Window: Shows visually the desired components, plots and 

perspectives. | 

e Main Menu: Allows control of the modeling functions of ANSYS, to include 


the preprocessor, the solution processor and the two postprocessors. 
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The user ıs able to modify the position and size of any of these windows and 
toolbars, therefore they will not shown here. When commands are referred to in the 
remainder of this thesis, they will be given in both the GUI menu-based commands and 
the direct text commands. An example of this 1s: 

Main Menu>Preprocessor>Material Props, etc, Typed Command 
where the “>” symbol indicates the use of a menu nested inside the previous one. To 
execute a desired command, there are usually several options available. For simplicity, 
these GUI commands will be given in a manner that 1s logical for the creation of the 
active model. In this thesis, only the commands that affect the model are included. 
Commands that affect the graphical display are omitted for spatial considerations. The 
syntax and specifics or the ANSYS commands used are detailed in Reference 25. 

The ANSYS preprocessor permits the creation and manipulation of a finite 
element model. This model can then be acted upon with loads and restraints, and then 
meshed into a finite element grid. The solver sub-program allows the selection and 
specification of different solution types, for example static analyses and modal analyses. 
For this thesis, a modal analysis was performed on the bare truss, and a transient analysis 
was performed on a controlled model. ANSYS has two different post-processors, one 
that works for single cases (or single steps of a transient analysis), and another that allows 
fusion of the entire transient picture. These en Е are called the POSTI and 
POST26 processors, respectively. 

When using the GUI for modeling and analysis, the user can also access a log file 
that gives the text commands that correspond to the graphic commands used. This can 
enable the user to recover a portion of lost data, or to easily recreate the model with 
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different parameters. The user can also access an error file directly from the GUI as an 


aid in resolving model problems. 


B. EXAMPLE MODEL 


In order to develop and demonstrate the ability to create an FEM of an actively 
controlled structure, and then be able to control it, a simplified model of a control strut 
was used. This simple model consisted of two piezoelectric elements, both of the same 
size and type, one to be used as a sensor, the other as an actuator. These elements were 
connected to a rigid base at one tip with BEAMA elements and subjected to an axial 
sinusoidal force at the other tip. The example FEM, called “piezo5” is shown in Figure 
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Figure 19. PiezoS Example Model (Screen Capture) 
This model was used only to verify that the algorithm worked properly for an 
actively controlled device, and was not intended to provide any concrete results. The 
sections that follow describe in detail how piezo5 was created and modeled in a transient 


environment, both uncontrolled and controlled. 
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1. Sensor and Actuator Integration 


The ANSYS element that has the required piezoelectric properties is the SOLIDS 
element. This element allows coupling of many of the physical properties permitted by 
ANSYS. For the piezoS demonstration model, two SOLIDS elements were required, one 
as sensor and the other as the actuator. To simplify the implementation, the same 
material and geometric properties were used for both elements. 

The SOLIDS element has eight nodes; therefore, linear shape functions were used 
to derive the local elemental matrices for use in the global FEM. Each of the eight nodes 
has up to six degrees of freedom: three-axis displacement, voltage, temperature and 
magnetic displacement [Ref. 23, p. 4-41]. It 1s of note that rotations are not included in 
the list of degrees of freedom for the SOLIDS element. 

The elements were integrated into the model by creating them adjacent to the 
structural elements that would support them. A rıgıd region was then created that would 
couple the degrees of freedom of the piezoelectric element to the end of the BEAM4 
element that formed the structure. The command line for a ngid region 1s: 

Main Menu>Preprocessor>Coupling/Ceqn>Rigid Region CERIG 

The rigid region was created only in the displacement degrees of freedom, as the 
SOLIDS element does not have rotations allowed at the nodes. Since the piezoelectric 
sensor and actuator are designed to detect and apply axial loads only, this was determined 
to be a sufficient modeling technique. A piezoelectric device is usually only shrouded, 
not constrained, in the off-axial direction. The rigidization of the ends of the SOLID5 
elements will modify its operating characteristics somewhat as the strain induced by 


fixing the cross-section of the element will reduce the strain applied in the axial direction 


46 


somewhat. Follow-on work should attempt to characterıze thıs modeling effect, or 
determine a resolution. 

Piezoelectric effects are, by convention, referred to the poling axis, which is 
usually the third axis of a Cartesian coordinate system. For the creation of this model, 
this convention was maintained. ANSYS allows the definition of a local coordinate 
system that will eier the specific alıgnment of an element. The local coordinate 
system can be created by reference to model keypoints, nodes, or the work plane origin 
[Ref. 22, CS]. The location is irrelevant so long as the alignment is correct. The work 
plane was relocated and aligned to the axial direction of the active element. The 
command line to create a local coordinate system aligned with the work plane origin is: 


Utility Menu>WorkPlane>Local Coordinate Systems> | 
Create Local CS>At WP Ongin = О РЕ 


To create the piezoelectric effect, a voltage difference needs to be applied to the 
elements. In order to do this a fixed voltage degree of freedom (DOF) was placed along 
one of the planes perpendicular to the poling axis. To preserve the ial symmetry, and 
to be consistent with the actual device, the nodes on the other perpendicular plane were 
constrained to be equal to one another by the command: 

Main Menu>Preprocessor>Coupling/Ceqn>Constraint Eqn er 

The final step in creating the model is the application of the specific piezoelectric 
properties. The SOLIDS element obtains the piezoelectric properties by the use of 
embedded tables. The piezoelectric table is a six by three table that represents the tensor 
between the forces and applied voltages, the [e] matrix. The table must be activated 


before setting the constants, and is referenced row by row, in a sequential manner (the 
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first row is 1-3, the second 4-6, etc.). Of note is that if any element has a voltage DOF, 
all or none must have a piezoelectric table. A suggested solution is to have very small 
piezoelectric constants for the non-piezoelectric elements. [Ref. 23, SOLIDS] 

The piezoelectric properties used for the PiezoS model were obtained from a 
verification model that was included in the ANSYS software package. The specific 
verification model was VM175, contained іп the directory ANSYSROOT\DATA\VERIF. 
VM175 contains the commands in text format for creating PZT-5A, an effective 
piezoelectric material. To prevent the model from searching for the magnetic and diei 
properties, in an attempt to solve these DOF, the relative permeability and thermal 


conductivity were set at 1.0. These added constants showed no effect on the results. 


2. ANSYS Parametric Design Language (APDL) 


It was desired to have a closed program, or a model that could perform an actively 
controlled data run, requiring no input from the user save to establish initial parameters 
and start the simulation. Two options for modeling the control law were explored. They 
are described in the paragraphs below. 

The first ле involved the use of COMBIN37 Elements operating on model 
node points to provide a “black box” addition to the model to simulate the controller. 
The COMBIN37 element has the capability of acting (expanding or contracting) on a 
derivative, double derivative or single integral of a node point value [Ref. 23, 
COMBIN37]. The element would be external to the model, fixed at one end, with a 
constant force applied. The stiffness of the element would be varied accordıng to the 


integral of the sensor voltage, given by a changing, potentially non-linear parameter. The 
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changing parameter is described by the term RVMOD,; its effect for integrating is given 
by the following equations: 


Bx M (4.1) 
РЕ. [Аха = [А sin(@t)dt = —K ro | (4.2) 
@ 
0 0 


С, 


RVMOD = RVAL +C,|CPAR|® +C,|CPAR 





(4.3) 
where RVAL is the current real constant value and RVMOD is the new value, when 
modified by equation (4.3). This leads to the output length of the element, shown here as: 

Ax =F /RVMOD. | (4.4) 

The changing length of the element would supply the input for a LINK11 element 
imbedded in the structure, serving as the actuator. The advantage for this model would 
be the speed of implementation, as there would be no external interfaces required by the 
user, save to modify the gain parameters. Due to the observance of the absolute value 
terms, this method abandoned as not immediately feasible for sinusoidal sensor input due 
to the mathematical implementation problem. Later work to allow operation on a bias 
would remove this problem for a sinusoidal signal. 

The second option, and the one that was eventually used, was ANSYS Parametric 
Design Language (APDL). APDL is a macro language that enables a user to run 
command streams from an external file. An APDL file is a text file that contains the 
commands as they would be typed into the ANSYS command interface with user 
determined variables [Ref. 22, /INPUT, *SET]. APDL was used in conjunction with the 
transient analysis capability to perform the control law application. ANSYS gives the 


user the ability to examine each step in a transient analysis; at the time of the step, the 
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APDL program can read the voltages from a sensor, perform a numenical control law 
function and apply a voltage to the actuator. 

APDL was originally created to allow the variation of model parameters in an 
optimization design. A capability is present to do automatic optimization, which was not 
explored at this time. APDL also contains some basic program control flow logic, such 
as 1f-then statements ña loops. These program control functions were exploited in the 
creation of a working controller. In order to use an APDL file in a loaded model, the 
command is as follows: 

Utility Menu>File>Read Input From: (select file), ¡INPÚUIA 

Up to 20 APDL subroutines may run, nested within one another [Ref. 22, 
[INPUT]. This number is reduced as loops are mere as ANSYS a level of 
nesting for them. This nesting utility allows for more complex — from within 
ANSYS to suit the user's desires. In this case, only one input file was required. Local 
parameters to these nested levels are possible by clearing the new variables before exiting 
the nested subroutine. 

To run the model, a routine was created that would automatically time-step 
through a transient analysis. A sinusoidal force was applied to the tip. This varying force 
was reapplied at each step with the use of an array-type variable, indexed to the time step 
number. 

Variables are created and deleted with the use of the *SET command. Vanables 
may be created as scalars, arrays and multi-dimensional tables. The sinusoidal applied 
force was created with an array that was sized based upon the length of time for the 


model run. The force was created as an array and as a one-dimensional table. The array 
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allows accessing by automatic indices. The tables are not pre-indexed, and cannot be 
accessed by the APDL program without programming the indices. However, for 
accessing the data in an output file, the tabular format provides simplicity; for this reason, 


both were used. 


3. Control Law Application 


The control law was implemented with the use of finite difference equations. In 
order to approximate the integrating function, a trapezoidal rule was used. The 


trapezoidal integral approximation is given by the following equation: 


iv = d £x 


t 


> (4.5) 


+1 


To perform this algorithm in the APDL macro, storage locations needed to be established 
for the single and double integrals, one for each. 

In simple runs of the model, it was found that the minute positive error from the 
double precision variables used by ANSYS caused a constant DC bias to be present, 
which destabilized the model. ANSYS tends to view the real number 0.0 as a random 
number on the order of 10°'. In order to prevent this occurrence, a simple digital high- 
pass filter was inserted in the system after the first integral term. This digital filter used a 
single pole and zero at 1.0 + 0.01 and 0.95 + 0.01 respectively. The transfer function for 


this filter is given by the following equation: 


Lee) 
a (20595) 7095) (4.6) 


Equation (4.6) is converted into a finite difference equation by the following 
equation [Ref. 25, p. 97]: 


S 


H(z)= —= 40 (4.7) 


Combining the previous two equations gives the finite difference equation to be 
inserted in the control signal path after the first integrating term, as shown in the 
following equation: 

Y(n)=X(n)-X(n-1)+0.95Y(n-]) (4.8) 
The frequency response of the digital filter is shown by Figure 20, which was generated 
with the MATLAB FREQZ command. As can be seen, there is not an appreciable phase 
error above 5.0- percent of the Nyquist frequency that is inserted by the use of the high- 
pass digital filter. The magnitude reduction is negligible above 2.0-percent of fe 
sampling frequency. The need for a filter in Piezo5 is consistent with the implementation 


of the dSPACE digital signal processor described in the previous chapter. 
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Figure 20. Simple Digital High-Pass Filter Frequency Response 


Finally, the control law can be obtained as a combination of the filtered integral 
term, and a double integral term, dertved from the first integral. These terms are 
multiplied by three gaın constants to get the control sıgnal to be applied to the actuator. 

The control law ıs given by: 


CS = (G,FINT +G,DBLINT)*G, (4.9) 


where CS, FINT_and DBLINT are the control sıgnal, the filtered integral and double 
integral terms, respectively. The term G, was used as a global system gain. 

After each time step is solved for, the single time post-processor was entered to 
obtain the voltages from the Sensor piezoelectric element ІТ the *GET command. The 
control law finite difference algorithm was performed and then applied as a voltage to the 
actuator element in the next time step. The control signals and the tip displacement were 
captured as well, for later analysis. 

The variable definition, load establishment, transient loop and control law 
integration are all included in the APDL macro “Truss.inp.” The “.inp” suffix was used 
to indicate that the file is used as an input file, and carries no special meaning with 


ANSYS. “Truss.inp” is included as Appendix G. 


4. Actively Controlled Model 

With the pieces assembled, namely the model and APDL macro program, 
“Truss.inp,” controlled runs to evaluate effectiveness of the control parameters could be 
performed. Initial trial runs of the model revealed that the initial start of the applied force 


created a transient response that operated at the natural frequency of the model. To 
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remove this effect and allow for examination of control of the driven frequency, damping 
was added to the model. 

ANSYS has three different damping options available, constant, Rayleigh and 
modal. Not all options are available with the different analysis types. For the transient 
analysis, only Rayleigh damping is available. Rayleigh damping creates a [C] matrix by 
using constants multiplied by the Mass and Stiffness matrices. The damping matrix is 
given by the following equation: 

[C]=0[M]+ £[K] (4.10) 
Rayleigh damping can be used to approximate the structure in the frequency region of 
interest, given the damping ratios for the given modes. The damping ratio, ©, is a 


function of the Rayleigh damping coefficients based upon the following equation: 


a fo 
cr p LAS 4.11 


A set of simultaneous equations for two frequencies can be set up, and then the 
Rayleigh constants may be solved. For this analysis, the damping coefficients set as 10- 
percent, giving « to be 0.085 and p to be 6.1e-4. The 10-percent value was arbitrary and 
was selected to damp out the transient vibrations. 

The time step was selected using the guidance presented in section 5.12 of 
Reference 26, the ANSYS Structural Analysis Guide. Under this guidance, a time step 
that was one-twentieth the applied frequency was used. This provided results that were 
accurate enough for this qualitative analysis. As discussed earlier, there is not an 
appreciable error above five-percent of the sampling rate. Sampling at one-twentieth the 


applied frequency will place the driving frequency at five-percent of the sampling rate. 
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The Piezo5 model was run under uncontrolled conditions and with varying G, 
settings. The model was excited for 2.0 seconds at 16.75 Hz, the frequency used, by 
Vlattas and Johnson [Ref. 14]. The 2.0-second run time was found to minimize the 
processing time while giving ample data for frequency analysıs. After initial trial and 
error, it was discovered that a polarity mismatch existed in the actuator. To correct this, 
gains of negative value were used. This solved the problem and enabled for reducing the 
amplitude of the tip motion. It was found to be more efficient to clear the database and 
start fresh prior to a new run to prevent corrupt constants from affecting the results. 

Nine data runs were performed on the Piezo3 model. The first was uncontrolled, 
and the rest varied the G3 system gain term from 0.0 to —20.0. The results are presented 


below. 






Run | G3 | Amplitude | .dB. — 
[1 | 0 | 1.670E-07 | 0.000 
| 6 | -8 | 1.530E-07 | -0.761 | 
| 8 | -13 | 1.497E-07 | -0.950_ 
ore ЛОС FO S 


Table 4. Piezo5 Active Control Results 
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Figure 21. Piezo5 Active Control Results 


(1 103X=%] 


VAL UE 





ROW 
Figure 22. Tip motion of Piezo5 


Figure 22 shows the motion of the tip of the Piezo5 model when the active control 
15 activated. The x-axis of the plot reflects the time index, each number represents a 
sample step; the active control is turned on at step 100. The amplitude is small, on the 
order of 10". As can be seen in the above results, the best case for the actively controlled 
model occurred when the G3 term was equal to -11.0 with a reduction in magnitude of — 


1.014 dB. The amplitudes are small, and the results are modest, but this 1s consistent 


56 


with the fact that PiezoS was intended as a demonstration model to prove the 


controllability of the FEM approach. 


С. NPS SPACE TRUSS 


1. Model Construction 


The first step in the ANSYS analysis ofthe NPS was to create a bare truss FEM. 
A simple structure, such as the NPS space truss can be rapidly constructed by using some 
of the features inherent to ANSYS. The WorkPlane is a user modifiable plane that И о 
the user to select points and draw on the model on the plane surface. In this manner, 
when a point ıs selected on the screen it ıs located where that point intersects the 
WorkPlane ın 3-D. The WorkPlane grid was a up such that the spacing between grid 
points was the distance between nodes on (s actual truss. Keypoints were located at 
these points by snapping them to the grid and joined by lines in same manner as the truss. 
The elements used to mesh the model were BEAM4s and MASS21s. Тһе ВЕАМ4 
elements used the geometry of the aluminum truss elements. The MASS21 elements 
were used to account for the additional mass of the fittings, connections, node balls and 


the knurled knobs that were attached to the node balls. 


2. Bare Truss Modes and Natural Frequencies 

To verify the FEM created with ANSYS, a modal analysis was performed. When 
an analysis is performed on ANSYS, the program will deliver a summary of mass 
properties and CPU allocation parameters. The mass of the model was compared with 


the mass tabulated in Appendix A. The mass of the FEM of the truss was 11.689 kg, as 
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compared to 11.708 kg for the actual truss. ANSYS was directed to provide the first 10 
natural frequencies and the first four mode shapes. 

The mode shapes obtained from ANSYS for the bare truss model are presented in 
Table 5. The frequency sequence follows the order and sequence of those obtained from 
the MATLAB FEM presented in Chapter JJ, and the experimental results obtained 
previously [Ref. 14]. 

ANSYS will provide to the user, if requested, various plots and tables of the 
results. The mode shapes for the first four modes were en in this manner. It was 
found that these mode shapes were identical to those obtained with the MATLAB FEM 


earlier. 
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Table 5. Natural Frequencies of ANSYS Bare Truss FEM 





D. INTEGRATION OF THE ACTIVE CONTROL ELEMENT 


The steps necessary to create an active finite element in a preexisting structure 
follow a similar pattern to the steps to create the sample model, Piezo5. In a realistic 
model, geometry and material constants become important. These step are laid out as 
follows: 

e Unmesh the FEM in the location er the active element. 

e Create geometry for the active elements: Lines for the connecting parts and 


volumes for the sensor and actuator. 
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e Create the material and mechanical properties for the connecting parts and the 
piezoelectric materials. 
e Orient the WorkPlane to align the 3-axis with the axial direction of the active 
elements and create a local coordinate system fixed to the WorkPlane origin. 
> Create Mesh: BEAM4 elements for the support connections, SOLIDS 
elements for the piezoelectric elements. Mesh volumes as a single SOLIDS 
element. Ensure that SOLIDS elements have the local coordinate system as 
the elemental coordinate system. 
e Create ngid region (in displacement DOF) to couple the ends of the sensor 
and actuator to the connecting pieces. 
e Apply fixed voltage DOF at one end of the sensor and actuator elements. 
e Couple the other end of the active elements in the voltage DOF. 
e Create the SOLIDS data tables for the piezoelectric properties (PIEZ). 
The NPS space truss active FEM has two different piezoelectric materials; the 
force transducer is quartz and the actuator is made of PZT. The constants needed for this 
analysis are the axial stiffnesses, piezoelectric constants and dielectric constants. The 


governing equations that ANSYS uses for the piezoelectric effects are in the stress form: 


15) b d "mis (4.12) 


For piezoelectric properties, ANSYS required the [e] matrix as opposed to the [d] 
matrix. By convention, the piezoelectric coefficients are given to populate the [d] matrix, 


and must therefore be converted to the [e] matrix by the following equation: 


[e] = [c] [d]. (4.13) 

Since the piezoelectric materials are only constrained in the axial (or 3) direction, 
some simplifications to the constant determination process may be used. Only the 33 
terms of the [d] and [e] matrices are important. By omitting the off-poling axis terms 
from the PIEZ data tables, the piezoelectric effect is therefore nullified, and the constraint 
that was presented earlier with the rigidization of the ends is alleviated. Therefore, 
Equation (4.13) can be simplified to 

€, 2 Ed. (4.14) 

For the quartz force transducer, the material specifications provided the sensitivity 
and stiffness (see Appendix A for the full component specifications). By using the 
stiffness and sensitivity, the d,, and E,, terms may be determined by direct combination. 
The elastic modulus was derived from the stiffness using the following equation: 

IOAN Je (4.15) 
The size used for the quartz element was a cube, 1.0 cm per side. 

The d33 constant for the actuator was obtained from the experimental work 
performed by Johnson and Vlattas in Reference 14. As the macroscopic properties were 
determined experimentally, it was possible to use a single element to model the actuator, 
even though it was composed of several PZT wafers, electrically connected. 

A summary of the material properties used for the SOLIDS piezoelectric elements 


is presented in Table 6. SOLIDS Material Properties. 
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Table 6. SOLIDS Material Properties 





A visual depiction of the piezoelectric sensor and actuator is shown in Figure 23. 


This figure was obtained from a hard copy screen capture from ANSYS. 





Figure 23. Installation of Active Members in NPS Space Truss 


- The LPACT was modeled by creating a MASS21 element with the mass of the 
device at the appropriate location. The transfer function was used to determine. the 
applied force. As was shown in Chapter II, the LPACT has a natural frequency at 


approximately 8-10 Hz. This location is less that the desired frequency range for the 
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truss, but as can be seen from Figure 11, there is large variation in the transfer function 
from 10 to 20 Hz. Using picked data points and the MATLAB function “Polyfit.m” a 
quadratic transfer function was derived. This function is given by the following equation: 


AMP = 0.2037 * FREQ? —7.0719* FREQ + 68.3564 (Ib/amp) (4.16) 


Using unit conversions and the electrical characteristics for the LPACT, a 
function could be applied in the macro APDL program that would simulate the change in 
the LPACT force that would D. applied in the FEM. The element that contained the 
LPACT was divided into two separate elements to create the location for the MASS21 
element and the location for the force application. 

The commands that were used to create the finite element in the NPS space truss 
were recovered from the log file. An edited version of these commands is included as 
Appendix H. Commands that controlled the display were among those omitted from the 
appendix. 

For the NPS truss active FEM, Rayleigh damping was again used. | Given the 
limitations of modeling Rayleigh damping for a structure with a large number of DOF, it 
was decided to select a value of 10-percent for the damping ratio, at the first two natural 
frequencies. Using the first two natural frequencies for the actual truss with the 
components installed, 13.113 and 17.125Hz, the Rayleigh constants, a and ß, were 
determined to be 9.1 and 1.1x107 respectively. The 10-percent value was an average of 
experimentally determined values in the band of interest [Ref. 14, p. 111]. Equation 


(4.11) was used in a series of simultaneous equations to perform this calculation. 


The APDL macro program, “Truss.inp” was modified to become "Act truss.inp", 
which served as the root program for the actıvely controlled NPS space truss model. The 
modifications include changing the node locations to the appropriate nodes in the truss 
model, and adapting the routine for the LPACT forcing function. The 20-to-1 gain 
provided by the Trek 50/750 Voltage amplifier was lumped into the SG term. 
“Act truss.inp” is included as Appendix I. 

The first test runs of the active NPS space truss FEM revealed a significant error 
in ANSYS. Once the model was complete, including the PIEZ tables, a modal analysis 
was attempted to determine the new natural frequencies, for comparison with the 
experimental data. ANSYS version 5.5 does not permit a modal analysis to be run with 
the PIEZ tables in the model. This path was abandoned and the transient analysis was 
attempted. Upon the attempt of the transient analysis an error was produced that stated: 

“Piezoelectric materials may not be used with the subspace eigensolver. Please 
use the reduced solver (MODOPT,REDUC command) with master DOF’s only at 
the displacement nodes." 

The subspace solver 1s an algorithm used in modal analyses to speed up the 
calculations, and was used for the both the bare truss model and the attempt at the active 
model. When the transient analysis mode is selected, ANSYS does not even permit the 
de-selection of the subspace solver. The problem was eventually resolved by exploring 
the model’s data base file using the command 

Utility Menu>File>Write DB Log File LGWRITE 

The LGWRITE command revealed that an artifact of the modal analysis run on 

the bare truss model was still present in the data base file. This was resolved by editing 
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the data base file to remove the references to previous modal analyses. This error has 
been brought to the attention of ANSYS, Inc and hopefully will be resolved in future 
versions. [Ref. 27] 

For the larger NPS space truss model, selection of a time step seemed to be more 
crucial due to the larger number of system degrees of freedom. A poorly selected time 
step can cause some numerical instability in a transient analysis [Ref. 15, p. 285]. 
Include the running of the model with negligible force. Early runs of the repaired model 
proved to be unstable after the third or fourth time step. By changing the sampling 
frequency, this problem was unable to be alleviated. Only by running with infinitesimal 
amplitudes could the model be run. 

This was later determined to be from —— data ın one of the data ӘНЕ 
Initially, the ANEL data table was used to define the anisotropic material properties. This 
table had incomplete data and was causing large shear stresses to exist in the model 
during transient analysis. By removing the ANEL data table, and placing the elastic 
constants іп an orthotropic state, this problem was alleviated, allowing the FEM to run 
through its specified time. [Ref. 28] 

The model was reconstructed from a version that was incomplete and had never 
had an analysis performed on it. After the bare model was reconstructed, a transient 
analysis was performed on the bare model using a version of the active APDL macro with 
the active parts removed by commenting. The active components were installed with the 
exception of the PIEZ data tables. The transient analysis was run successfully on the 


model. The PIEZ data table and material properties were created and the model was able 
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to run without instability. Further evaluation of the NPS space truss FEM is contained in 


the next chapter. ^ 


pe ALTERNATE SIMULINK MODEL 


In order to evaluate an active piezoelectric control element in a Simulink 
analytical model, its operating characteristics must be modeled. The derivation of this 
modeled transfer function, from force sensor to piezoelectric actuator is complex, and 
will be explained in detail. 

For this model, built from the independent FEM used in Chapter II, it was 
assumed that the Een element and the strut are the same material and length and 
co-located. This 1s not realistic for a quantitative model, but it is sufficient to prove 
controllability of this model. The embedded force sensor detects motion from the axial . 


expansion or contraction of the element. The nodal displacements must be converted to 


the local frame by 
{x} =[DCM] {x}. (4.17) 
For a truss element, only the axial change in length 1s important. This change is given by: 


бы=и-щ=[-1 0 .. 1 0 5). (4.18) 


Using the stress-strain relation gives 


F Ж (4.19) 


axial T Le 
The active element can be modeled as springs in series; the effective stiffness of the 


element is therefore the reciprocal of the sum of the reciprocals, given by the following: 
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From Newton's first law, the force in the sensor is the same as felt by the member. 
Therefore, the change in length for the sensor 1s given by 


F 
RECTA (4.21) 


Sensor 
sensor 


The sıgnal voltage output from the sensor is given by the following equation 
o (4.22) 


Sensor Sensor 


The control law used, as described previously, 15 presented below 


C = +) V one l 
5 5 


(4.23) 


where C is the control signal K, is a constant term reflecting the gain from an amplifier 
on the output of the force sensor to the processing computer. 

The output from the controller is applied through a voltage gain amplifier and is 
applied through the piezoelectric effect to the long axis of the element as a control force. 


This relation is given by 


-1 

0 

| : 
Қаны) = 1 ey (4.24) 
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where K, is the amplification gain from a voltage amplifier and the computer output gain. 
The voltage limiter in the dSPACE model was not included in this model to simplify ıts 


construction. 
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Converting to the global coordinate frame from the local frame gives 


ЈЕ 


control 


 =[РСМ] {Ешь}. (4.25) 


Combining the above equations yields 


-1 


0 
: G G k iezo а; 1020 
apo = [ремт“ қ К, ae ШЕ 1240 ][DCM]" {x}° (4.26) 
5 sensor  33Sensor 
| 0 | 
| | 
Equation (4.26) ıs cumbersome; as most terms are constants, this equation can be 
reduced to 


¡E 


control 


| =[АРЈК (S45 иво · (427) 


Where [ISO] is the vector that isolates the desired displacements from the truss converts 
them to the local frame and determines the change in element length. [AP] is the 
application vector that takes the axial output of the active strut and applies it to the truss 
structure. For this analysis, the exact value of the constants is unimportant, and served 
only as an alternative exercise for the author. Therefore, all the constant terms were 
lumped together as the variable K. This allows the computer to perform the integral and 
double integral transfer functions on a single scalar, as opposed to the entire state vector. 
The active control model was built in Simulink, using 300 for G, and 100 for G,. 
This was in keeping with the results from, Reference 14. A twenty-second run was 
performed on each location, acting on an impulse disturbance acting at node three in the 


y, z and torsional x-directions. This three-degree of freedom impulse was selected to 
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excite several different vibrational modes of the truss. For consistency, for each run, the 
output was obtained from the y-displacement of node 3, and a lumped gain, K, of 30,000. 

A time step of 0.0001 sec was used, as larger time values were found to be 
numerically unstable in the model. Kwon [Ref. 15, p. 285] states that the time step must 
be less than the minimum model period divided by x. An analysis of a desired time step 
based upon the maximum frequency of interest criteria was not performed as stability was 
the only desired objective. Any attempt at qualitative analysis should include an 
appropriate time-step, based upon ES frequency of other appropriate criteria. 

Early evaluation revealed a bias in the values for sensed displacement, similar to 
the experimental results obtained earlier. This is believed to be an artifact of the 
numerical analysis, and was alleviated by inserting a Butterworth high-pass filter in the 
control stream with a high-pass frequency of 3 radians/second. The complete Simulink | 
model for this analysis 1s shown in Figure 24. 

In order to prepare the FEM for the Simulink edel the program *NPS prep.m," 
included as Appendix J, was created. "NPS prep.m" served to perform the matrix 
inversions and multiplications required and to create the [AP] and [ISO] matrices. 

The model produced a qualitative validation of its controllability. In the test runs, 
the control law was activated at 5.0 seconds by changing the lumped gain value from zero 


to 30,000. A screen capture of the sensed control input 1s shown in Figure 25. 
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Figure 24. Simulink Active Controlled Model 


ad Control Еогсез 





Figure 25. Sensed Signal From Active Simulink FEM 


69 





N COMPARISON OF ANSYS FINITE ELEMENT MODEL TO NAVAL 
POSTGRADUATE SCHOOL SPACE TRUSS 


A. METHODOLOGY 


The purpose of the experiments performed in this thesis was to verify the active 
FEM developed in the previous chapter. To do this, five experiments were performed on 
the NPS space truss to evaluate the IDIFF controller’s performance. 

Before performing control parameter experiments, the truss was excited’ with 
random noise from a digital spectrum analyzer to verify the natural frequencies. The first 
experiment that was performed was an evaluation of the controllability of the first natural 
frequency. The second natural frequency was targeted for the next experiment. For the 
third Sen a repeat of the experiment performed by Johnson and Vlattas [Ref. 14] 
was done. The fourth experiment was to combine the first two natural frequencies with 
variable phasing. Finally, random noise was applied and controlled with the IDIFF 
controller. 

For each experiment, a series of runs were performed with varying gain settings. 
The value of the read sensor was used to help find the optimum point. The goal in these 
experiments was to reduce the vibration on node 26, simulating the location of a sensitive 
piece of equipment. The intent of the LPACT was to provide excitation from a simulated 


vibrating component operating at critical frequencies. 


pel 


B. EXPERIMENTAL CONTROL SYSTEM 


1. Experimental Setup 


In order to detect the motion of the truss, four Kistler 8690C10 three-axis 
accelerometers were used. Two of the accelerometers were place at the extremities of the 
truss, nodes 26 and 41. FEM modal analysis revealed that these locations experienced the 
largest motion when excited at the natural frequencies. These tip sensors also will 
capture the motion from the first several modes. The second pair of accelerometers were 
placed at the midpoint of the truss arms, on opposing sides to the tip accelerometers 
(nodes 49 and 18) to capture the higher modes where the ends of the truss do not vibrate 
as much as the middle. 

The output of the accelerometers was fed into a Kistler Piezotron Signal 
Conditioner, model 5124A. The signal conditioner does not атты the signal, and only 
served to power the accelerometers and filter the output. The accelerometers experienced 
significant dnft during the course of the experiments. This drift was able to be easily 
isolated from the signal of interest with a spectrum analysis. 

The - of the digital control system was the DS1103 digital signal processor 
from dSPACE. The dSPACE system received inputs from the quartz force transducer, 
after conditioning to perform its control function. The conditioned accelerometer signals 
were fed into the ADC inputs for structural monitoring. 

To perform a pseudo-integrated monitoring function a Hewlett-Packard HP 
54601А four-channel oscilloscope was used to monitor the read control sıgnal from the 


force transducer and the output signal to the piezoelectric actuator. The oscilloscope was 
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set up such that both channels could be monitored at a sufficient resolution to obtain a 
rapid qualitative feel for the behavior of the system at a glance. 

To supply the excitation to the LPACT, a Hewlett Packard HP 33120 digital 
signal generator as used. The HP 33120 has the capability to provide many different 
wave shapes, at a variety of frequencies and amplitudes. Selecting a frequency and 
amplitude can be Ет” с by direct keypad insertion at the front of the signal analyzer, 
or by using a dial to scroll through the selected significant digit on the display. 

All the connections between components were provided y varying length BNC 
connector cables. These cables had no special specifications. The completed 
experimental control system can be seen in Figure 26. Expenmental Layout, shown 


below. A detailed description of the cable connections 15 provided in Table 7. 
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Figure 26. Experimental Layout 
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Table 7. Experimental Cable Connections 





Once the cables are installed, the truss may be activated. The first step in 
activating the experiment is to power up all the equipment that will be used. Finally, the 
Newport Vibration Isolation table will need to be floated by opening the valve to the 
Nitrogen es Over time, the compressed gas will bleed from the pneumatic cylinders 
supporting the table; therefore, care must be taken to ensure that the valve is shut upon 
completion of the experimental data taking, to prevent the exhaustion of the Nitrogen 
supply. As the experiment is located in the basement of Halligan Hall, this evolution 
takes some time to perform, and draws away from the time available for experimentation. 

The dSPACE ControlDesk program is then called up on the host PC, using the 
steps described in Chapter III. То summarize, the ControlDesk experiment 
“5Т controller" is opened, which will open the panel layout file and the trace variable 
file. Finally, the *ST controller 3" program is loaded into the CPU. 

At this time the signal generator may be activated and tuned to the desired 
frequency and amplitude. The LPACT does not activate until the LPACT switch is 
placed to “enable.” When enabling the LPACT, caution must be taken to ensure that the 


actuator is being held secured, this will prevent a large transient from going through the 
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truss and damaging sensors due to a large response. Also, when active, adjust the applied 
voltage slowly using the vernier dial on the front of the signal generator, to prevent the 
same transient. The frequency may be adjusted by the use of the keypad while the 


LPACT is enabled, as this does not cause a noticeable transient to occur. 


2. Modal Verification 


Before the experiments were performed to evaluate the active control system with 
the IDIFF controller, the modal frequencies were verified with a Hewlett-Packard 
HP35665A Dynamic Signal Analyzer (DSA). The HP35665A has to ability to generate a 
random noise signal within a desired frequency range and can analytically isolate the 
desired frequency range. This analysis was also performed to determine the highest . 
effective frequency, which was used to obtain the time step size for the Simulink model. 

The DSA was set up in two-channel mode, with a split screen to allow viewing of 
both channels. Each channel was set up for fast Fourier Transform (FFT) analysis. One 
channel was set to record magnitude; the other channel was set for phase measurement. 
Тһе modes become apparent where there are local peaks in magnitude with a 
corresponding 180-degree shift in the phase. As random noise was used to obtain the 
frequency response, averaging was used to smooth the signal. The averaging method 
used for this analysis was vector averaging. 

Vector averaging preserves the complex values of the signal, and allows 
averaging phased data, such as an FFT analysis [Ref. 27, p. 4-313]. The phase 
information in the vicinity of the modes was still degraded due to the rapid phase shift at 


these locations. By looking at the data in mid-average, it was easier to determine where 
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the modes were. The HP 35665A DSA setup procedure is included as Appendix K. HP 
35665A DYNAMIC SIGNAL ANALYZER SETUP. The equipment setup is shown in 
Figure 27. The random noise source, in addition to being the driver for the truss, was fed 
back into channel one of the DSA as transfer function input. The response from the force 


transducer was sent to channel two as transfer function output. 
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Figure 27. HP 35665A Modal Determination Setup 

The DSA ıs also capable of saving the data on a 3-1/2 inch floppy disk. These 
data files are saved ın a format called SDF (Standard Data Format). A routine was 
provided with the HP 35665A DSA to perform a conversion to standard MATLAB 
MAT-files. This routine is called SDFTOML. This routine is executed from an MS- 
DOS wındow with the following syntax: 

MS-DOS Prompt> SDFTOML <SDF filename.EXT> <MAT-filename.MAT> 

The data that is recorded in the MAT-file consists of four variables, 0211, 0211x0, 
0211x1, and 0211xl. These variables are the data, which is a matrix of dimension two by 
the number of resolution lines, plus one, and the frequency specifications. The variables 
0211x0, 0211x1 and о211х1 are the starting frequency, the frequency resolution and the 


ending frequency, respectively. 
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The captured data was transferred to the dSPACE host PC ın the sub-directory: 
C:\Space_Truss\truss99\experiments\HP DSA. Six runs were performed, examining the 
frequency at varying resolutions over zero to 200 Hz. By examining smaller frequency 
ranges, a greater resolution was achievable. The data files were examined with the M-file 
DSA PLOT.M, included as Appendix L. The range from zero to 50 Hz is presented as 


Figure 28. 
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Figure 28. Frequency Response from Random Noise Response (0-50Hz) 


The first and second natural frequencies will be targeted for this thesis, as well as 
the 90-degree offset point examined by LT Vlattas and LT Johnson [Ref. 12]. Due to the 
resolution of the DSA, the precise location of the natural frequency was not obtained, but 
a small band was used with the vernier feature of the signal generator to lock in on the 
natural frequencies. 

As can be seen from the results there is a new natural frequency in the 8-10 Hz 
range. This frequency is due to the addition of the LPACT to the truss structure. This 
natural frequency is local to the LPACT strut, given the low stiffness of the LPACT drive 
spring, and was therefore not included as part of the frequency range of interest. This 
was confirmed by examining the first mode shape from the truss with the active 


components installed, as generated in ANSYS. 


py 


The frequency range of zero to 200 Hz was examined in order to determine the 
maximum frequency of interest for the truss. Ав сап be seen from Figure 29, the 
magnitude of the frequency response begins to taper off at 200 Hz. For this reason, 
frequencies above 250 Hz, or sampling at 500 Hz, should provide for negligible aliasing 


effects. 
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Figure 29. Frequency Response from Random Noise Response (0-200Hz) 


C. EXPERIMENTAL PARAMETERS AND RESULTS 


1. Data Capture 

The dSPACE software package also includes two function libraries that allow 
real-time access to the DS1103’s memory. These two utilities are called MLIB and 
MTRACE. These functions were designed to be executed from within the MATLAB 
workspace, or from m-files that the MATLAB program is running. MLIB provides 
access to the parameters that the CPU is using for the running dSPACE application, 
MTRACE provides for continuous or timed data capture capability [Ref. 30, p. 6]. A 


pictorial description of this relationship is shown in Figure 30. 
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Figure 30. MLIB/MTRACE Relationship [From Ref. 30, p. 6] 


The MLIB and MTRACE routines are invoked by function calls from within 
MATLAB m-files. The specific function desired 1s called as a parameter to the function 
call, with the details of the command being secondary parameters. Examples of this are: 


mlib(‘DesiredCommand’,’parameter name','parameter value','etc."); 
mtre1103(‘DesiredCommand’,’parameter_name’,’parameter_value’,’etc.’); 


The MLIB and MTRACE libraries were used to instruct the DS1103 CPU to 
provide 20 seconds of data for each experimental run, at the Simulink model-sampling 
rate of 500 Hz, for a total of 1,000 samples. The parameters that were obtained were the 
accelerometer ae: the force transducer input and the dSPACE signal output to the 
piezoelectric actuator. The program “Acq data.M,” included as Appendix M, contains 
the code that was used. “Acq data" was modified for each experiment to ensure that the 
data obtained would be deposited ın the correct subdirectory. 

Using the guidance contained in E 14 and 30, the program was 
instructed to capture and process the data in the following manner: 


e The CPU was polled to venfy that “ST Controller 3.ppc" was running. 
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e Variables were defined and dSPACE instructed to find their memory 
locations. - 

e Capture settings were established, to include sample length. Data can either 
be stored on the DS1103 local memory or sent to the host PC for storage on 
the hard drive. 

° After the data capture was complete, it was sent to the MATLAB workspace 
as a matrix of dimension, number of variables by number of samples. 

e The data were E into individual vectors for plotting and performing 
frequency analysrs. 

One of the improvements in the ControlDesk program was the inclusion of a 

D&D function that can retrieve the variable name directly into the MATLAB editor for 
inclusion 1n an M-file. When a Simulink model 1s created and sent to the CPU with the 
RTI program, it creates a trace file that contains the memory locations and names of the 
model parameters. These parameters are sorted by dSPACE by function and block 
location, as described in Chapter III. The D&D function saves time by alleviating the 
need for the programmer to re-type the entire variable name, as seen by dSPACE, and 
reduces programming errors. 

After the first run for each experiment, which was uncontrolled, the program was 

altered to use the first set of data as a comparison. All the captured data was saved as a 
separate MAT-file for later processing. The data processing routine is included as 


"Data proc.m," Appendix N. 
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2. Experiment 1, First Natural Frequency, 12.50 Hz 


The first experiment performed was a control of the first fundamental frequency 
of 12.50 Hz. This frequency was verified prior to experimental runs by using the vernier 
dial on the HP 33120A Signal Generator. The vernier dial was adjusted until the largest 
response was seen on the oscilloscope display. The amplitude on the signal generator 
was varied until a noticeable and clearly sinusoidal response was seen on the 
oscilloscope. An amplitude of at least 1.0 Volts peak-to-peak (Vpp) was required to 
obtain this response. | 

Raising the voltage to 1.2 Vpp allowed the laser diode to provide a clear visually 
qualitative picture of the controllability of the truss. The laser diode reflected the side to 
side nature of the mode shape by its projection Bh the laboratory wall. 

The first set of runs was еа phases. First,.the single integral term, 
IGain, was set at 300 and the system gain, SG, was modified to find the limit of stability. 
Once the limit was reached, SG was backed off slightly to allow for variation in IIGain, 
the double integral term, to find an optimum point. This was then repeated with IIGain 
being the dominant variable. Finally, the amplitude was lowered to attempt to determine 
if the 100 mVpp signal reported in Reference 14 was possible. 

For each run, the amplitude and power reduction from the uncontrolled truss at 
node 26 in the Y-direction were determined. Node 26 is at the opposite extreme from the 
truss as the LPACT. These values are listed for all the data runs in Appendix O. 

In order to conserve space, the graphical plots of the captured information (control 
signal input, dSPACE control output and the three-axis response of the four 


accelerometers), presented in Appendix O, will represent a snapshot of the runs from all 
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experiments, showing a generic uncontrolled case, unstable cases that are borderline 
stable and rapidly unstable and the best case controlled case. These plots are presented as 
Figure 46 through Figure 56. 

During the test runs, it was discovered that the Nitrogen had run out of the bottle, 
requiring a replacement. The next few data runs were slightly inconsistent with the first 
set. Due to the somewhat random nature of the data, there is some scatter from one 
measurement to the next, and from day to day. In order to obtain a consistent series of 
results, an experiment should be completed as soon as possible, within the same day, or 
hour as the data runs are only 20-seconds long. 

When instability occurred, the active strut traveled from saturation on one side to 
saturation on the other at a frequency that is different that any of the truss’ natural 
frequencies. High frequency transients are observed as the actuator enters and exits the 
saturated region. The spectral results show this transient as a series of spikes on the PSD 
plot. Also observed was accelerometer drift that occurred during the experimental runs. 
This was a low frequency effect and was apparent on the spectral analysis, but appeared 
in the lower frequency band (<2 Hz), and therefore was not of further interest. 

In all, 28 trials were performed, in addition to the uncontrolled case. The best 
results were obtained from trial 15, with parameters of IGain-300, IIGain--400 and 
SG=1.75. The reduction in magnitude from this trial was 18.54 dB. 

The relative contribution of the IIGain parameter was two orders of magnitude 
less than the single integral term. This relationship can be simplified by the following 


equation: 


А. 
| А сов па = — эт cot (Sa) 
0) 


It was discovered that the double integral term, in order to be amplified to a level 
where its effect could be felt on the truss, drove the system more rapidly into instability 
by amplifying the noise as well. When used in conjunction with the [Gain term, the 
IIGain term can be most effectively used to adjust the frequency response when damping 
out the disturbance. The IGain term refers to rate feedback, and affects the damping; the 
IIGain refers to position feedback, and modifies the system natural frequencies. 

Of note from the results 1s that as the desired location 1s reduced in amplitude, the 
other locations monitored show an increase. This is due to the installation of only one 
control strut in the truss. In order to target multiple degrees of freedom, and different 
mode shapes, more control struts would be required. This would also boost the control 
authority and further reduce the magnitude. Care must be take in selecting a second strut 
location to ensure that the desired effect ıs obtained. 

If multiple modes are targeted from multiple control struts, a method must be 
developed to assist with the selection of the active strut locations. One method is to use a 
weighing of the strain energies (in a manner similar to that in Chapter I, but for multiple 
modes) to select the optimum locations. This method was used successfully in previous 


research to control a cantilevered truss [Ref. 31, p. 394]. 


3. Experiment 2, Second Natural Frequency, 13.81 Hz 


The second set of experimental runs was targeted at the NPS space truss’ second 


natural frequency, 13.81 Hz. As with Experiment 1, this value was obtained by adjusting 


83 


the vernier dial on the HP 33120A Signal Generator. This experiment proceeded in a 
manner similar to experiment one, and comparable results were obtained with 21 trials: 
The best case achieved was in trial 8, with. parameters: IGain-300, IIGain —300 
and SG=1.75, at a reduction in amplitude of 19.02 Hz. This reduction was not at node 
26, as was used in experiment 1. For this mode the controller was unable to reduce the 
amplitude of node 26’s y-directed motion. Good results were obtained from node 41, the 
diagonal counterpart of node 26 at the other extreme of the truss. This trial was the best 
overall case during the course of research and testing. Graphical results from this trial are 


presented in the figures below. 
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Power Spectrum Magnitude (dB) 


Node 41 - X. Y Z-Axis - Time Data 
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Figure 33. Exp. 2 Trial 9 Node 26 Response and PSD 
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4. Experiment 3, Repeat of Vlattas And Johnson, 16.75 Hz 


The third experiment was a repetition of the test performed by LT Scott Johnson 
and LT John Vlattas [Ref. 14]. This experiment searched for a target frequency that had a 
90-degree phase offset from the actuator to the sensor. In their controller, band pass 
filters were used (high pass filters are currently used) to target a specific frequency. In 
this case 16.75 Hz was found to have a 90-degree offset. Use of the DSA to examine this 
region with greater frequency resolution confirmed that this was the location of the 90- 
degree offset. 

With the exception of the filters and the new dSPACE system, the experimental 
setup was identical. However, when starting the system, it was found that a 100 Vpp 
signal, as reported in Reference 14, did not produce a noticeable signal from the force 
transducer. Once the amplitude was raised to 1.0 Vpp a measurable signal appeared. 
Further increasing the amplitude of the generated disturbance signal to 1.2 Vpp allowed 
the use of the laser diode in qualitative study. 

Sixteen trials were conducted in this experiment. The best results that were 
obtained were during tnal 8, with gain parameters of IGain-300, IlGain=-100, SG=1.75 
giving a reduction of 17.93 dB. These results were consistent with the approximately 15 
dB reduction achieved previously [Ref. 14, p. 98]. A different value of IIGain was found 
to have the best response during this experiment than the runs performed by Lieutenants 
Johnson and Vlattas. This 1s possibly due to a slight shift in the frequency response 
owing to the moving of the truss from the Newport vibration isolation table during the 


salvage efforts after the flood in early 1999. 
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5. Experiment 4, Two Modes With Variable Phasing 


In this experiment, the first two fundamental natural frequencies, 12.50 and 13.81 
Hz were targeted EO to observe the broadband performance of the controller 
and to verify that the superposition principle is in effect. According to the superposition 
principle, linear signals may be separated, have linear operations performed on them and 
then be recombined to produce a single output. 

To create the dual frequencies, the dSPACE controller, created 1n Simulink, was 
modified to include a dual-frequency signal generator, with the capability of varying the 
starting phase of one of the signals. This phase variation was intended to explore the 
superposition idea just discussed. The ControlDesk experimental layout was also 
modified to allow quasi real-time adjustment of these parameters. The modification to 


the Simulink model is shown in Figure 34. 
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Figure 34. Dual-Frequency Modification to ST_Controller 
For this experiment, twelve trials were performed, at six different phasings. In 
each case, an uncontrolled run was performed, and a controlled run, operating with the 
parameters: IGain-300. IIGain-100 and SG-1.75. The phase offset for the second 
frequency was varied from zero to 240-degrees in 60-degree increments. The results 
from the six cases were all similar, and resulted in a reduction in amplitude of the highest 


peak of 7.29 to 9.26 dB. This reduction was in the first frequency, 12.50 Hz. The second 
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frequency, 13.81 Hz showed a reduction of about 15 dB. Each case showed similar 
results, verifying the superposition principle as applied here. The beating phenomenon 


from two close sinusoidal signals, and their spectral response are shown in Figure 35 and 


Figure 36. 
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Figure 35. Exp. 4 Trial 51 Node 41 Response 
An entire representative case of the dual-frequency experiment is shown in 


Appendix O, Figure 57 through Figure 59. 


6. Experiment 5, Random Noise 


The final experiment performed was intended to further assess the broadband 


performance of the truss controller with random noise. The truss was excited with white 
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Power Spectral Density - Controlled vs Uncontrolled (Tnal 51) 
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Figure 36. Exp. 4 Trial 52 Power Spectral Density 

noise, generated from the HP 35665A Dynamic Signal Analyzer. The DSA white noise 
function has a band-pass filter inherent in the output that makes the majority of the power 
in a desired spectral region. The region selected from this experiment was 10-35 Hz, as 
this contains a large number of modes. Outside of this band, the noise input tapers off 
linearly with respect to the logarithm of the magnitude [Ref. 28, p. 4-213]. The random 
signal generated from the DSA was echoed to the oscilloscope. Examination of this 
signal revealed that it was not truly random, but consisted of a series of sinusoidal 
fragments of randomly varying frequency. 

Twenty trials were performed. The first ten trials were performed uncontrolled, 
the second set of ten trials were performed with the controller active, with the parameters: 


IGain=300, IGain=100 and SG=1.75. The sets of runs were processed with the data 
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processing routines and their PSD’s averaged to find a mean value of the controlled and 
uncontrolled case. Examination of the data revealed a variation of +/-5dB. 

The change of the highest peak of the uncontrolled cases to the highest of the 
controlled cases was 6.69 dB. This number is not truly representative of the broadband 
performance of the system. The PSD is shown in Figure 37. It shows that the response is 
significant in the targeted range, and reductions of 10-25 dB were observed in the 
frequencies targeted for the first three experiments. There 1s an increase of about 10dB in 
the natural frequency that was added by the installation of the LPACT (about 8 Hz). This 
shows that the addition of the LPACT, especially the integrated spring, is significant to 


the frequency response of the truss, and should be examined in further experiments. 
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D. ANSYS FEM 


1. Model Setup and Initiation 


As discussed in Chapter IV, the SOLIDS elements are not permitted to have a 
modal analysis performed on them if they use PIEZ data tables. Also, if a modal analysis 
were performed on a model before the inclusion of the PIEZ tables, an error in the 
database would ensue that would prevent the running of subsequent transient analyses. 
To work around this issue, a copy of the active NPS space truss active model was used to 
delete the PIEZ data tables and obtain the modified truss natural frequencies. This model 
was not used for any transient analysis. Future versions of ANSYS will permit the PIEZ 
data tables to be used in modal analysis [Ref. 27]. It ıs not known ıf the transient 
problem will be resolved. 

The natural frequencies obtained by this analysis were significantly different from 
the original modal frequencies. This is 1n part due to the added mass of the ГРАСТ апа 
control apparatus (The mass 1s approximately 20-percent above the bare truss mass). The 
natural frequencies obtained from the modified FEM are presented in Table 8. 

Frequenc Frequenc 
9.187 13.113 
11.299 17.125 
11.743 29.688 
17.781 72.075 
26.332 89.371 
6 | 33.010 96.153 
59.797 103.046 
i8. a n EA 110.338 


| 9 | 77.830 135.820 
81.508 140.448 


Table 8. Natural Frequencies of Truss With Active Components 












E 


Upon examining the mode shapes from the postprocessor, it was determined that 
the first natural frequency was located at the LPACT strut and was due to the 2.5 kg mass 
vibrating there. This mode was not present in this form in the real truss as the LPACT 
has a stiffer adapter to the truss than the relatively long and flimsy aluminum beams in 
the FEM. However, in the real truss, the LPACT has a spring that does have a natural 
frequency of about 8-10 Hz (as seen in Appendix A). This frequency was seen in the 
frequency determination with the DSA and subsequent random vibration experiment, 
previously in this chapter. 

The second mode, at 11.299 Hz is a combination of the effects of the LPACT 
addition and the original first mode. Examination of the mode shapes revealed that the 
second mode combined the properties of the LPACT with the first original mode. The 
third modified mode shape was observed to be a combination of the first and second bare 
truss mode shapes. This produces a diagonal motion of the truss, and is consistent with 
the experimental results qualitatively observed with the laser mounted above node 52. 

When comparing the ANSYS modified natural frequencies with those 
experimentally obtained from the work of LT Johnson and LT Vlattas [Ref. 14, p. 110] it 
can be seen that there is a large discrepancy in the frequencies. It is believed that the 
previous experiment did not target the natural frequency that was added by the LPACT, 
and that the two analytical modes at 11.30 and 11.75 Hz were molded together in the 
experimental analysis. 

As stated in Chapter IV, the SOLIDS element is connected to the structure with 
the use of rigid regions that constrain the element from contracting as it expands in the 
axial direction. This will induce an artificial stiffness into the model if not examined. 
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The NPS space truss FEM was modified to delete all of the components not directly 
connected to the active strut, leaving only the SOLIDS elements and a few BEAM4s. А 
static analysis was performed, setting a 100V signal on the end of the piezoelectric 
actuator. 

As reported ın previous work [Ref. 14, p. 25-28] and shown in Figure 7, the 
actuator has an expansion of 50um for an applied voltage of 100V. Under the static 
analysis, the actuator had a displacement of only 37.1-um, 74-percent of the expected 
value. The inverse of the 74-percent results shows that an increase of 35-percent is 
required to match the experimental data, roughly the value of Poisson’s ratio. The value 
of d,, in the PIEZ data table was increased to 9568 from 7100 to account for this effect. 
This patch was only possible because of the prior work that measured the exact operating 
properties of the piezoelectric actuator. 

The actual truss active element was preloaded with shims to give a compressive 
bias for the force transducer and piezoelectric actuator. In the ANSYS model, this effect 
was not included due to the mathematical modeling of the SOLIDS elements. A limiter 
was included, but was set to +/- 60V. 

The runs were designed to run for 4.0 seconds of simulated time. As the time-step 
interval of 1/20" the driving frequency was used, a limit was discovered that affected the 
final time used. ANSYS will only permit the performance of 1000 transient steps in an 
analysis. This required the reduction of the total time for the 16.75 Hz data set. If longer 
data runs are required, ıt would be necessary to either save the final conditions and apply 


them at the start ofa new data run or use a larger value for the time step. 
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The active control APDL macro, “Act _Truss.inmp,” was modified to allow the 
expanded use of AÁPDL”s capabilities. - The macro was divided into two programs, one 
that set up the initial parameters and saved the calculated data, and the other to perform 
the transient analysis. In this manner the workstation was directed to perform 10-40 data 
runs at a time (limited by the time desired for completion and analysis by the author). As 
each data run was approximately 15 minutes and 40 seconds, this allowed use of the 
computer processor in the off-hours. 

The files saved by the overarching APDL macro were ASCII files that contained 
the node 26 y-direction motion, sensor input and controller output. A MATLAB 


algorithm was developed to process the data in a rapid manner. This algorithm 15 


included as Appendix P. PROC_ANSYS.M. 


2. ANSYS Series 1: 11.30 Hz 


In this series of data runs, the second natural frequency, 11.299 Hz was targeted. 
As discussed above, this frequency is the combination of LPACT effects and the first 
original natural frequency. The overarching APDL macros were directed to examine 
variations in IGain, [[Gain, and finally to attempt to get a better performance by ШЕ... 
slightly off the IGain limit and attempting to find a best operating condition. Detailed 
electrical schematics for the sensor and actuator were not available, so trial and error to 
determine the polarity of the control loop was performed. It was observed that negative 
values for SG provided control, indicating a polarity correction ın the sensor of actuator. 

The ANSYS model displays different instability characteristics than those of the 


NPS space truss. The voltage limiter drives the control signal into saturation, which 
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causes transients to be seen as disturbances on the amplitude graphs and as side lobes on 
the PSD chart. Depending on whether the model was limited by IGain or IIGain, 
different results were observed. The IIGain term causes larger side lobes on the order of 
80dB. Examples of these effects are shown in Appendix Q. ANSYS RESULTS, Figure 
61 and Figure 62. The IIGain term causes a different settling time constant that the IGain 
dominated cases. An example of a slightly saturated data run is presented as Figure 63. 

In all, 27 data runs were performed for the first series. The best overall results 
were obtained with trial 16, with a reduction in node 26 amplitude of 22.4521 dB. This 
value was the best obtained for the five different series of data runs. This data run is 
presented as Figure 38 and Figure 39. Trial 16 was run during the examination of the 


double integral effects alone, and represents a position only style of feedback. 
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Figure 38. Series 1 Trial 16 Results 
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Power Spectral Density - Controlled vs Uncontrolled - Tnal 16) 
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Figure 39. Series 1 Trial 16 Power Spectral Density 
From the above figure, an increase can be observed at the 33-Hz frequency range, 
at the location of the sixth analytical mode. As can be observed, there is an increase 
outside the range of second and third natural frequencies, the ones targeted for specific 
analysis. This effect may be due to the Rayleigh damping term used that is effective in 


the desired range, and gives poor broad frequency range damping characteristics. 


3. ANSYS Series 2: 11.75 Hz 


The second series targeted 11.75, the third natural frequency. Due to the decaying 
nature of the LPACT signal as frequency increases, smaller amplitudes were observed for 
the uncontrolled case in this series, as compared to the first. Consequently, a slightly 
larger gain was applied to the controller to reach saturation. At this frequency, the model 


displayed similar saturation transient characteristics as the series performed at 11.30 Hz. 
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Ihe APDL macro ran the same series of gain constants as the first series. 
Additional runs were performed at the end to find a limit with both gain terms. Of the 34 
runs performed, the best results obtained were in trial 34, showing a reduction of 18.0605 
B. This i was at the saturation limit. An unusually low value was observed for trial 
22, reflecting that the gain constants were in error. The APDL Macro has already been 
modified when this discovered, but the consistency in the nearby tnals supports this 
conclusion. Care is required when programming the APDL macro. Future studies could 


examine the use of easy to read data tables to minimize errors. 


4. ANSYS Series 3: 16.75 Hz 


Series three was intended to observe the experiment performed by LT Johnson 
and LT Vlattas [Ref. 14]. Initial runs on this series revealed that very different results 
would be obtained. As the negative gain was applied, an increase in the amplitude was 
observed. An example of this was trial 3, shown in Appendix Q as Figure 65. Further, 
the controller drove itself into saturation at gain levels that were far below those of the 
first two series (the first series became saturated at 2150 and —125 values of SG, senes 
three was saturated at 50). As the now positive SG was raised, the model entered a 
saturation region and did not re emerge to go to the other side, as was the case ın the 
previous two series. 

Adjusting the three gain settings did little to correct this occurrence. The limit 
between stable and unstable was difficult to find. The best results obtained were in trial 
4, showing a reduction of 19.6464 dB at the target frequency. This run, however, caused 


the amplitude of the signal in the 8-12 Hz band to increase by 100 dB, a significant effect. 
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Trial 4 1s shown as Figure 66. Trial 5, the saturation and instability limit, is shown in 


Figure 67. 


5. ANSYS Series 4: Dual Frequencies 


The dual frequency experiment was repeated in this series of data runs. The 11.30 
Hz and 11.75 Hz frequencies were used. As these two frequencies are close in value, a 
long time beating phenomena was expected. For this reason, and due to the good results 
obtained with series one and two, the sampling time was raised to Y 10% the highest 
frequency, with a sample run of 8.0 seconds to allow the beating phenomena to be 
observed. 

These trials were run with gain settings that were satisfactory for both frequencies 
to prevent instability. As the only noise in a digital system is due to round off error, these 
results were much more consistent than the parallel experiment performed on the actual 
NPS space truss. For all trials, results of about 15 dB were obtained. A typical 


controlled trial is presented in the appendix as Figure 68. 


6. ANSYS Series 5: Random Noise 


For the final data series, twenty runs were performed, ten uncontrolled and ten 
controlled. The sampling frequency was set at 200 Hz, at 4.8 seconds to maximize the 
collection within the 1000 time step limit imposed by ANSYS. ANSYS has a random 
number generator function in APDL that was used to obtain the random signal. 

The PSDs obtained were averaged, and their values compared. The reduction 


from the peak uncontrolled to the peak controlled value was 6.0371 dB. The PSD is 
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presented as Figure 40. It can be observed that the originally targeted frequencies, in the 


11-13 Hz band show reductions of about 15 dB, which is consistent with series four. 


Power Spectral Density - Controlled vs Uncontrolled - Trial 20) 


-195 


BO P. Neo t 


№ 
=> 
O 


' 
N 
N 
о 


Power Spectrum Magnitude (dB) 
N 
ол 


-225 pun 


-230 ee fA 5. 





5 10 15 20 25 30 
Frequency (Hz) 


Figure 40. Senes 5: Random Power Spectral Density 


The purely random noise in this series did not produce as large a reduction that 
the first two series due the filtering action. The digital filters required three steps of 
consistent data to get a good data trend. The purely random nature of the disturbance 
signal did not give ample time for the filters to perform properly. Both the ANSYS 
results and those obtained from the actual truss show an increase in the 8-10Hz range 


owing to the LPACT contribution. 
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E. COMPARISON AND SUMMARY OF RESULTS 


The best cases from each run, experimental and analytical are presented here, as 


Table 9. 










Experimental | Analytical_ | IGain | liGain | SG_ | Results (dB) | 
| Exp.1Trial15 | - | 300 | -400 | 1.75 | 18.5410 — 
НЕ ОКЗ l o шн а 
__хр.2 тла | № * | 300) 300”/| 1157)” 19020248 
____-- _____ Бепе5 2 гап 34 | 300 | 40,000 | -100 | 18.0605 | 
ПОЗЕ Е | = | 5007 о еш 
J | Series3rund | 300 | 0 | 25 | 196464 — 
__Бхр. 4 Тпај51 |___-- | 300 | 100 | 1.75 | 92564 | 
| |Series4rn31| 300 | 0 | -75 | 15.5659 | 
Г вр б -- |300) 100 ООО 
E ATA Ses O o h a СО И 


Table 9. Summary of Results 
















The results obtained with the analytical model in ANSYS behaved in a manner 
similar to those obtained by experiment. Overall, the amplitudes in the analytical model 
were slightly higher than those in the experiment, by about 2-3 dB. The system gain 
term, SG, was inclusive of the 20-to-1 Trek 50/750 voltage amplifier that was present in 
the actual space truss. This reflects most of the difference in the values of SG when 
comparing the analytical results with the experiment. When this factored in, the 
analytical model has the ability increase SG by a factor of three before saturation. 

When the tests were driven into instability, a difference was discovered between 
the ANSYS model and the actual truss. The actual truss would become very unstable, 
and le at the saturation limit for a period that was much longer than it was driven. 
The analytical model would ~ enter saturation, and change states at each period, at 
the same frequency as the driving signal. Both showed high frequency transient 


responses at the saturation lımits. 
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The actual truss and the ANSYS model exhibited an overall reduction ın the 
frequency range of interest, and an increase in the magnitude at other frequencies, 
especially at the LPACT frequency range (8-10 Hz). Overall, the ANSYS model, while 
more stable owing to its saturation vice instability showed poorer response in the higher 
frequency bands. A more detailed model, including better damping effects, may better 
capture this effect. 

Finally, the analytical results were more consistent overall than the experimental 
results. This was due to the computational nature of the experiment, and could not 


include the random effects that any real system would experience (i.e. noise). 
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УІ. CONCLUSIONS AND RECOMMENDATIONS 


=ч 


By using IDIFF control for the Naval Postgraduate School Space Truss, vibration 
reductions of 15-20 dB were realized with the use of the dSPACE data acquisition and 
processing system. The ANSYS actively controlled finite element model showed similar, 
qualitative results, amplitude reductions of 18-22 dB. This demonstrated that the method 
devised to actively control an FEM is valid, and can predict the control authority in a 
structure with active components installed. 

The ANSYS model was devoid of system noise, except that which occurs with 
small round off errors that ле present ın computer processors. This may be the cause of 
the slightly better results and different stability and saturation characteristics observed in 
the experiments and data runs. 

The analytical method developed using ANSYS, which is commercially available, 
has room for refinement. The modeling of the actuator with a SOLIDS element using 
rigid regions to attach it to the structure caused some inconsistencies in the initial data 
test runs. A calculational patch was used to remedy this problem, but further analysis is 
required to refine the method for integrating an active control de into a ruce, 
finite element model. 

The responsiveness of the NPS space truss to the IDIFF controller was positive 
(1.e. reduced the amplitude) in the examined degrees of freedom. Other DOFs showed 
large increases in amplitude. The addition of more actuators to the NPS space truss is 
recommended to allow better vibration reduction in the truss as a whole. Other 


controllers could be explored that would lead to larger reductions. An option for a 


103 


prospective controller is to use the MATLAB finite element model contained herein to 
create a dynamic observer using the accelerometer input to create a full state feedback 
controller for the truss. 

A feedback-style controller is only effective to a limit (as the sensor signal lowers 
it reduces the controller output) and 15 not usually able to fully eliminate the vibration in a 
structure 1f continuously driven. A feed-forward controller, such as the clear box system 
developed for the Ultra-Quiet Platform, another active vibration control experiment in the 
SRDC, may be implemented as a means to accomplish this. These controllers rely on 
some knowledge of the disturbance, either by sensing directly or by estimation, and are 
limited in that regard. 

The purely random nature of the signals created in the ANSYS data series 5 did 
not allow a good look at the spectrum, as the random nature prevented the filters from 
being able to have the proper effect. In order to demonstrate a good capability at all 
frequencies, a frequency sweep algorithm should be written with APDL to perform an 
analysis over the range of the frequencies desired. This is a task best suited for APDL. 
As each run was about 15 minutes in length, an entire data set would take days to 
perform. 

The method created within of using ANSYS to perform an active control 
simulation on a structural FEM can be used on any manner of structure. It is able of 
being added to any pre-existing ANSYS model, or created with a new model. In this 
manner large vibration-susceptible structures, such as the International Space Station can 
be examined with active vibration reducing devices to determine their effect on the 


structure when in orbit. 
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APPENDIX A. NPS SPACE TRUSS CHARACTERISTICS AND COMPONENTS 


Mass Properties ofthe Bare and Modified Truss: 


Bare Truss | Mod Truss | Masses (kg) | Truss (kg) | Truss (ka 
TRAS 0 ia CO OO 2.276" | 
Act. Strut#1] 0 =| 1 | 0290 | 0000 | 0.290 | 


Total Mass 11.708 14.118 
Truss | 


Table 10. NPS Space Truss Mass Properties [After Ref. 14] 













Bare Truss Natural Frequencies 


Mode Number | Modal Testing |MATLAB FEM | ANSYS 
Ref. 14 


[3 p M54 [448 | 3428 — 
[| æn | 3287 | 28895 


oe | 745 | 7224 | 7298 _ 
[8 | 30101 | м | 9661 — 
[— 9 | 623 | 2021 | 11541 - 


Table 11. NPS Space Truss Bare Natural Frequencies 









Table 12. Truss Element Properties 
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COMPONENTS 


Piezoelectric Translation Model P-848-30 S/N 


Volts Microns Microns Microns Percent 
O oom ошо Т 20007 ГЕ 
5.98 
10.12 
zo 17206) ОО шш 
БЕСІ” ПЕТ = у A 
_ 90 | 4405 | 4607 | 202 | 41 __ 
| 1000 | 4914 | 4914 | 000 | 0,00 __ 


Table 13. Expansion and Contraction Data for Model P-843.30 [From Ref. 14] 

















Open Loop Travel (0-100V) 45 um +/- 20% 
Closed Loop Travel | 45 шп 
Stiffness 33 N/um +/- 20% 
Force Generation (Blocked) 1500 N +/-20% 
Push/Pull Force Capability 800/300 N 
Torque Limit (at tip) 350 mNm 
Capacitance 5.4uF +/- 20% 
Dynamic Operating Current Coefficient 15 LA/(Hz-um) 
Unloaded Resonant Frequency 10 kHz +/- 20% 
Operating Temperature Range -20 to +80 °C 
Mass (w/o cables) 536 
Length 73 mm 
Table 14. P-843.30 Operating Characteristics [Ref. 32] 
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Planning Systems Incorporated LPACT [Ref. 18] S/N CML-030-020-1 


Force Constant (K,) 5.5 lb/amp 


Max. Current 1 атр 


~3 % (or critical) without force loop, 
up to =a. u with force loop on 


(Stoke ww A E 2 | OQ inches По 


force 


Gravity Offset Spring Rate 
clamp and accelerometers) 
LPACT Model (low frequency) Output Force(Ib) 
(refer to Figure 2 for measured uc ANN 2. 


FRF from current to force of 
LPACT) #@ 5” 


s+2£0,5+ 0, 


Servo Amp Model aie Cumentlemp) _ O oj р 
Servo Command (V) Я | 


Force Loop Model Servo Amp Voltage Command (volts) 
(see section 3.3 for definition of Pr oof Mass Accel (g) 


terms) 2 
De Ка K force > 


Current Command(amp) 


2 
(S + Wore )(S + Wy) 


Rate Loop Model Servo Amp Voltage Command (volts) 
(see section 3.3 for definition of Primary Accel(g) 
terms) d 5 


KKK Wiese 


rate rate” r 


M V. LG 4 wu) 


Table 15. LPACT Characteristics 
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Cable Assembly cable | connect to LPACT Electronics connect to LPACT 
(all on rear panel) Component 
Black Coax "To Coil 6" Blue Pigtail from coil 
(banana plug to BNC adapter) (BNC) 


Blue Coax marked with ‘From Secondary Secondary Accelerometer 


Red Tape Accelerometer’ (BNC) on Proof Mass (microdot) 
Blue Coax ‘From Primary Accelerometer’ Primary Accelerometer 


(BNC) on Co-Locate Ring (microdot) 
Table 16. LPACT Electronics Connectivity Guidelines [From Ref. 18] 
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Figure 41. LPACT Control Electronics Rear Panel [From Ref. 18] 
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where m = mass of proof mass (2.9 Ib) 


Figure 42. LPACT System Level Block Diagram [From Ref. 18] 
PCB Piezotronics Model 208B02 General Purpose ICP Force Sensor S/N 15021 


Sensitivity 50 mV/lb (11240 mV/kN) 
(Specification) 

Sensitivity (Measured) 50.80 mV/lb [Ref. 33] 
Dynamic Range -10 lb to 100-1b 


Stiffness 1.0 kN/um 
Temperature Range -54 to 121 *C 
Sensing Element 





Table 17. PCB Model 208B02 Operating Characteristics [Ref. 17] 
PCB® Piezotronics Type 484B Signal Conditioner S/N 2086 


Notes 


Unity Gain 
Set CPLG to DC & Bias to 6 V 


Kistler Instrument Corp. Accelerometers: 
(Note: g = 9.807 m/s”) 


Type Senal Number + х-ах15 + у-ах15 + Z-axis 

8690C10 C112398 495 490 494  mV/g 
8690C10 C112399 487 490 490 mV/g 
8690C10 C112400 499 500 494  mV/g 
8690С10 (112401 497 491 505 mV/g 


Kistler Instrument Corp. Signal Conditioners (Multi-Channel Couplers): 


Type Serial Number 
5124A (twelve channel) C74930 


О 


Trek Voltage Amplifier: 


Type Serial Number Notes 
Trek 50/750 none Required calıbration on 10 September 1998. 


Mini-cal to verify setting and linearity 
Performed by author on 2 July, 1999 

Two channels that can be used with both active 
struts. 


dSPACE GmbH, Germany, S/N 3192 


@ Motorola PowerPC running at 300 Hz 
Memory 128MB RAM 
Input Channels 20 ADC 


Ouput Channels 8 DAC 0-10V 


Hewlett Packard HP 33120A Signal Generator S/N 
settings 


Hewlett Packard HP 54601 A Digital Oscilloscope S/N 3134A02713 
Four Channels 
Frequency response up to 100 MHz 


Hewlett Packard HP 35665A Dynamic Signal Analyzer S/N 
hookups and settings 


SOFTWARE DOCUMENTATION 


[dsPACERIING —— 71393777 рт" ы 


ANSYS/Multiphysics 












Table 18. Software Documentation 


114 


APPENDIX B. FINAL ROOT NPS.M 


~ 


Final root.m modified from 

Example gg 

Modified by LT Carey M. Pantling for ME4613 Final project 
eigenanalysis of NPS Space Truss 

Bulids systems mass and stiffness matrices 


до до dO go до до 


Variable descriptions 
coord = global x,y and z coordiates of each node 
nd = nodal connection vector 
k = element stiffness matrix 
kk = system stiffness matrix 


до до до до до 


% m = element mass matrix 
% mm = system mass matrix 
% Ef = system force vector- 


index = a vector containing system dofs associated with each element 
bcdof = a vector containing dofs associated with boundary conditions 
bcval = a vector containing b c values associated with 

the dofsS in 5ecdor! | ` 


- = | on oT чж «ж oT oT oT «ә. «ғ» «ж» -т- а «-. «ж U «ше ан. ee ан. = = = = = = = = = = = = = = = = = = = = = 


до dP № до oe 


% 


clear; 

nel=161; $ number of elements 

nnel-2; %$ number of nodes per element 
паоЕ=в; $ number of dofs per node 
nnode=52; $ total number of nodes in system 


total system dofs 


до 


sdof=nnode*ndof; 


load nps coord dat и о ут 
load nps_node.dat % 161 by 2, 161 elements, nodel node2 
load nps_bcdof.dat % 24 by 2, 24 degrees constrained 
coord=nps_coord; nd=nps_node; 
bDedot-nps bedor, 1); bceval=nps_bedof(:,2); 
el=70E9; % GPa elastic modulus 
area=1.9686982E-5; $ m^2 cross-sectional area 
xpr-1922122175-107 $ m^4 moment of inertia of cross-section 
rho=2800; $ kg/m^3 mass density per volume 
ff-zeros(sdof,1); $ initialization of system force vector 
kk-zeros (sdof,sdof); $ initialization of system matrix 
mm-zeros (sdof,sdof); $ initialization of system matrix 
index=zeros (nel*ndof,1); % initialization of index vector 
for iel=i:nel % loop for the total number of elements 
index=feeldof (nd (iel,:),nnel,ndof); % extract system dofs 
% compute element stiffness and mass matrices 
{k,m] =feframe3 (el,xi,area, rho, coord (nd(iel,1),:),coord(nd(iel,2),:)); 


kk-feasmbli(kk,k,index); % assemble each element matrix into system 
mm=feasmbli(mm,m, index) ; 
end 


Jui 


% add the concentrated masses at the nodal translation dof 

for iel=1:nnode 
Start= (161-1) *ndor-+1; 
fini=start+2; 
+ ада lumped masses for the node points 
mm(start:£lni,start:fini)=mm(starc: fim, Stamens 0 0663) ~Cye erm 
$kg added to each node as concentrated mass 

end 


€ apply boundary conditions 
ГКК тт, ЕЕ] =Ғеартус Кк, ит, ЕЕ bcdof, bevel ) ; 


(V,D] =eig (kk, mm) ; 
рео  - слали; 
fsol=sqrt({fsol)/2/p1; 
% purge NaN's 
fsol(1snan(tsol))=1! 
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APPENDIX C. FEFRAME3.M 


function [k, mi=teframe3 (el,x1,area,rho,xyzl,xyz2) 


From feframe2.m by Young Kwon 
Modified for 3D by LT Carey M. Pantling 
Purpose: 
Stiffness and mass matrices for the 3-d frame element 
nodal dof (u 1v 1 w 1 theta 1 theta 2 theta 3 

u2 v 2 w 2 theta 4 theta 5 theat 6] 
element stiffness per FEM/MATLAB by Kwon, etc. p 264 
adds lumped masses at ends for the screws and fittings 


AP AP AP AP 99 99 99 99 99 


Synopsis: 
[k,m]=feframe3 (el,xi,leng,area,rho,xyzl,xyz2) 


до ор оо 


Variable Description: 
k - element stiffness matrix (global) 
m - element mass matrix (global) 
el - elastic modulus 
xi - second moment of inertia of cross-section 
area - area of beam cross-section 
rho - mass density (mass per unit volume) 
ху21 - coordinates of first node (1x3) 
xyz2 - coordinates of second node (1x3) 
kl - element stiffness matrix (size of 12x12) (local) 
ml - element mass matrix (size of 12x12) 
dx,dy,dz - differences in the three axes 
leng - element length 


99 99 99 99 99 99 99 92 90 92 99 909 99 


% ul,vl,wl - local unit vectors (global vectors i,j,k) 
% OA- off axis vector for local unit vector definition 
% a,b,c - constants for the stiffness matrix 


J=rotational inertia term =2*x1 
G=shear modulus = E/2.6 assumes nu=0.3 


до до 


% compute elemetal dimensional data 
AXexyz2 (I) хут (1), ЧЕМУ (2). 
leng=(dx*2+dy*2+dz*2)%0.5; 


compute element rotation matrix 

local aligned with x along axis, y,z orthogonal 

since tube segments, y,z do not matter 
DCM defined as (x)g=9[c]1*(x)1; g[c]l is the transformation matrix 
ul=[(dx/leng), (dy/leng), (dz/leng)]; 
cr=eross (ul 10,011). 
mag= (cr (1)%2+cr(2)%*2+cr(3)*%2)%0.5; 
if mag-=0 % if {ul}x{i}-=0 then 
OA=[0,1,0;-1,0,0;0,0,1]*ul'; 2 rotate z-90deqrees 
else 

OA=[0,0,-1;0,1,0;1,0,0]*ul'; % rotate y-90 degrees 


до оо 


до до 


епа 
ул=стовв (ПОД) у уй=у1/ (у1 (1) о ИО) ОЕЕО ВИО Е 
wl=crossiul ee wl=wl/ (will) ы (27721 (37720055: 


n7 


DEM= ["ulr(1)Z I ти ет ои 
ul(2) ; У м/н 
а (З) са ee eel 

r=zeros(12); 

r(123,1:3)=DEeM: 

r(4:6,4:6)=DCM; 

r(7:9,7:9) =DCM; 

r(10:12, 7@e72)=DCM:; 

Tri; 


% stiffness matrix at the local axis 
kl=zeros (12);; 
G=el/2.6; 

J=x1*2; 
al=el*area/leng; 
a2=J*G/leng; 
b1-12*el*xi/leng^3; 


axial compression 

axial torsion 

pure bending due to deflection-y 
b2=6*el*xi/leng*2; pure bending due to deflection-y and rotation z 
b3=2*el*xi/leng; bending due to rotation z 

cl=bl; $ tube therefore symmetric 

c2=b2; % tube - 

c3=b3; % tube x 
kl(1:6,1:6)-[a1, OF 0:0 (а-а 

"AE 0, [Do O "DE a EE 


оо d? o9 o? o? 


0 
OF О, СВ Um CIE EM e 
OSOS 0, ал OO wae 
0; 0- C27 mO а ве вене" 
e do О, Оо ТУУ Ва ЈЕ 
КІ(1:6,7:12)-|-аТт/ 0 О or = 
0,-Бі; ~i ОТО ОТ 
OF 0v c1. Ос о 8 
м ОУ са ағалы С. 2 
0% Оро 2“ Оэ а 
0, -b2, Ver) ae NISI „ 
К в D2 
I 7212-7312) = 057 "ово о ENSE 
DS DIE О О: 
o7 О СІ 0 c2 MENO 


9, 97 ас, 0, De. 
0, A CA e О а 
0, -bD2, 0, o Ө 2 БОЕ 
% stiffness matrix at the global axis 
Ја = екет 


t compound mass matrix per MSC/Nastran 
ml-zeros( 12,12), 

mass=rho*area*leng/12; 

mi(1:3,1:3)=5*ey8 13). uml 7:5 v9 ml Mc E 
ml (1:3,7:9)=eye(3); mI S amma 
ml=mass*ml; 

t€ add mass for the fittings, with no change in stiffness 
ml(1:3,1:3)=ml(1:3,1:3)+0.015*eye (3); 
m1l(7:9,7:9)=ml(7 93,7 2) ОВЕ БЕ ЕЕС 

% mass in the global system 

msr mI T: 


118 


APPENDIX D. NPS MODES.M 


ao 
oo 
до 


%%8%%%58%%%%%%%%%%%%%%%%%%% % % SESS 
NPS MODES.M Plot code for mode shape plotting % 
By LT Carey M. Pantling % 


A до oo o9 oe 
о С 
oo 0 
oo (D 
до N 
ao 
de =] 
ado ы. 
do t 
до f 
oe LJ 
99 | 
о қ 
до О 
до О 
oe ct 
ae | 
до B 
до "8 
A? Y 
ao ° 
% 3 
ae 
o? Fh 
99 О 
X K 
oe 
до Б: 
до (ї 
99 0 
oo 
% ©’ 
oo f» 
(D 
de 'O 


+ plot the modes 

Ва: 

атр=2.0; 

тойе-у(:,259); 

disp=nps_coord; 

for index=1:nnode 
start=(index-1)*6+1; 


disp(index,1)=disp (index,1) +amp*mode (start); % х-соога 
disp (index, 2) =disp (index,2)+amp*mode (start+1); % у- 
disp (index, 3)=disp (index,3)+amp*mode (start+2); % 2- 


end 


$ isolate the real components of the eigenvectors 
disp=real (disp) ; 


% form the longeron combinations 
plots (disp(1:2, 1) disp T? а а тор 2 о 
hold on 


piot3(nps coord(1:2,1)9BDnpsNesord(r]:2,3) 00D NcooOrd(1:2 2) E 


view(-30,20); 
ах15([-2 3 -1 2 -1 21); 
plot3(disp(3:14,1) disp(3:14,3) disp(3:114 Oo za 


Plor Empszeosrdt > РАИ) пре Собака ЛА З ips coord Nu E 


view(-30,20); 
axis([-2 3 -1 2 -1 2]); 
plot3 (disp (15: 261m ATsp(15S 26 есер ис 262) IE EN 


РОЗ (пр "coomd (15. 26,1) mps cora (MISS Mp Seo oa 26/2) 


view(-30,20) ; 
axis([-2 3 -1 2 -1 21); 
plot3.idi1sp(27:28,1) .dispi27 ВИЗА Бр ОВОО ИрИ 


plót3 (nps coord(27:28,1) Mps eoore 27, 258 o) прашев zu Zen, 


view(-30,20); 
axis([-2 3 -1 2 -1 2]); 
plot3 (disp (29:40,1) disp(29:140,3) disp(29:140,2),'b*="); 


plots (Nps e895rd1292.409 1), NpsreooranN 29.20, 3) ,0pszevord 297210725 


view(-30,20); 
acs а. 
pLOES (disp (41-5271) , disp (41-523) disp (41522) ee 


pLot3(nps coord (41-5251); Mpsuceord (4125275) пр-во. 


view (-30,20); 
axis([-2 3 -1 2 -1 2]); 


% Cross Cies 


rorsıindex- 12 26 
line Gia а bla ae 


MIO 


Paley, 


"qna 


a: 


UE 


К); 


11пе(2,:)-алвер іі пдехе267:); 
plot3 (line(:,1)% іле 24 Іле (:,2) шағаны; 
епа 


$ verticals 
torsındex= 1112227 28] 
line (7, 2)=disp (index, :); 
line (2, :)=disp (index+7,:); 
plot3(line(:,1),lyYBpe ame UENIT EE 
end 
for index= 3:14 
iinet >) =disp (index, =); 
linet: )=displindex+i2 7); 
purors (line(=,1), line(: 73), пе: ЭУ D-n); 
end 
for index= 29:40 
linei- disp(index i. 
11пе (2, 2) =4а15р(іпаехті2 5 
plot3 (line(*,1) lineas), lme (rob. 
end 
Поше нове 
title('Mode 4'); 
print -djpeg mode4ml.jpg 
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APPENDIX E. NPS SPACE TRUSS MODE SHAPES 
The first four mode shapes obtained with MATLAB are presented first. The deformed 
shapes are shown with the undeformed node locations for comparison. The first four 
mode shapes generated by ANSYS are presented second, with both the deformed and 


undeformed shapes. 





Mode 3 Mode 4 
Figure 43. NPS Space Truss Mode Shapes with MATLAB 
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The mode shapes that were generated from ANSYS were directed to show both 


the deformed and undeformed states. The undeformed states appear as dotted lines in the 


following figures. 
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Figure 44. NPS Space Truss Mode Shapes with ANSYS 
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APPENDIX F. NPS STRAIN.M 


$ Strain energy determination, by LT Carey M. Pantling 
Uses Final root nps for the determination of the system 
Variables required in workspace 


до до 


% kk - global stiffness matrix 
% cr - global force vector 
% mm - global mass matrix, not used except as part of program 


bcdof - DOF for boundary conditions 
beval - values for the BCs at DOF above 
coord - matrix (nnode by 3) of nodal coordinates 


99 оо oe 


% па - matrix (nel by 2) of nodal connections 
% е1 - Elastic Modulus 

% area - cross sectional area 

% Program Variables 

$ disp - displacements from solution 


delsq - elemental change in length, squared 


до 


$ get the strain energies for the mode eigenvectors 
disp=coord; % reset displaced positions 
amp=1.0; 
mode=V2 (:,260); 
% get the new nodal positions 
Тән index=-1:nnode 
start=(index-1)*6+1; 


disp (index, 1) =disp (index, 1)+amp*mode (start) ; % x-coord 
disp (index, 2) =disp(index,2)+amp*mode (start+1); % y- 
disp (index,3)=displindex,3)+amp*mode (start+2); % z- 

end 


% get the elemental strain energies U=0.5EAd*2/L 
for index=1:nel 
index; 
% original element length 
dx=coord (nd (index,2),1) - coord (nd (index,1),1); 
dv-eoeordind (index, 2), 2) -esoraind (index 1) 20 
dz=eserdind (index, 2), 3) -zsserdmdtindez, 3); 
leng= (dx*2+dy*2+dz*2)*0.5; € original length 
$ new length 
delx=disp (nd (index, 2) ,1) -disp (nd (index, 1),1); 
dely=disp (nd (index,2),2) -disp (nd (index,1),2); 
delz=disp (nd (index; 2) ,3) -dispimd (index, 1) 73). 
newleng= (delx”2+dely*2+delz*2)*0.5; % new length 
change=newleng-leng; 
str eng (index) =0.5*el*area/leng*change*2; 
end 
[str,in]=sort (str_eng) ; 
% flip the matrix 
for index=1:nel 
strout (index) =str (162-index) ; 
inout (index) =in(162-index); 
end 
Bead (streur (i300) ) | 
insut (1.30). 
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APPENDIX G. TRUSS.INP 


| ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЕЖЖЖЖЖЖЖЖЖЕЖЖЖЖЖЖЖЖЖЕЖ 


! * Truss Control APDL Program for ANSYS version 5.5 
| * For simple sinusoidal disturbance 

!* Wnitten by LT Carey M. Pantling 

! * With assistance by Sheldon Imaoka 

! * Last Modified 16 July 1999 


| ЖҰЖЖЖЖЖЖЖЖЖЖЖЖ Ж ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЕЖЖ Ж 


! First Load Truss model with mesh and BC's 


Y KK k k k k k K K K K K K K K K K K kK K K kk K ok k K k K K K K k Æ kK k 


! * Define Variables 
Yoko ak ok ok ok ade ale ok ale ade ade ole ok ok K K K K K K Fk de ae ad OK ade al ode ale ale ade ade ad OK OK OK OK OK OK OK OK OK ok OK OK Ok OK KK 


BSEDBISRO 16.75 ! disturbance frequency 

ЗЕКЕ оО ! control sıgnal to piezo, ınitial zero 
*set,PV ‚node(0.300,-0.004,0.004) ! Piezo voltage node for application 
*set,RS, ! read signal from piezo, defined only 
*set,S V,node(0.150,-0.004,0.004) ! Sensor voltage node for detection 
*set,G1,300 ! Control gain constants, see below 
*set,G2,100 

*set,G3,0.0 

*set,OLDRS,0.0 ! the old RS 

*set, OLDINT,0.0 ! previous integral, IC 0.0 

*set, OLDFINT,0.0 

*set, OLDDBL,0.0 ! prev DI, IC 0.0 

*set,S T[RTTIME,0.0 ! Start time = 0.0 seconds 

* set, FINITIME, 1.0 ! Finish Time in seconds 
TIMESTEP=1/FREQ/20 ! TIMESTEP 


*set,Pl,acos(-1) 
* set, DISPNODE,node(0.375,0.0,0.0)! DISPNODE = the node where 41 is in truss 


“set. CSCheck ! tracks the control signal 
*set, RSCheck ! tracks the read signal 
finish 


| ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ Ж 


! * Set up load disturbance, at least 20 steps per period _ 

| ЖЖЖ ЖЕТ А та жаа жат ЕТЕКТЕ ЕЗ жа 
*set, NUMSTEP,nint(FINITIME/TIMESTEP) 

*set,FORCFUN, 

*set, FFP, 

*dim,FORCFUN,array,NUMSTEP ! array for forcing function 
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*dim,FFP,tableNUMSTEP ! table for plotting 

*do,ICOUNT,1,NUMSTEP 
FORCFUN(ICOUNT)=sin(2*PI*FREQ*((ICOUNT-1)*TIMESTEP)) 
FFP(ICOUNT)=sın(2*PI*FREQ*((ICOUNT-1)*TIMESTEP)) 

*enddo 

*set, NODE41Y 

! creates an array to get the node 41-y displacements 

*dim, NODE41 Y,table. NUMSTEP 

*dim,CSCheck,table, NUMSTEP ! gets the control signal 

*dim,RSCheck, table, NUMSTEP ! copies the sensed signal 

*dim,INTcheck,table, NUMSTEP 

*dim,DBLcheck,table, NUMSTEP 

*dim,FINTchk, table. NUMSTEP 


| K K KK K kK kK k kK K K KKK K kK KK K K K K K K K KK K K K K K K K K K K K K K K K K K K K K K K K K K 


! * Do loop for loading at each time step 
[sb ok oe eoe me ook koe ook ok oe ve ok ok КЖ БЕ ЕЕЕ ЕЕ а а E E EE О СА E E E 
*do,ICOUNT,1,NUMSTEP ` 
! Set current time (actually, end of current step) 
*set, CURRTIME,STRTTIME+ICOUNT*TIMESTEP 


/solu ! enter the solution processor 
ЖГІСОСХМТ е ел 
antype,trans,new ! Start new transient analysis 
*else 
antype,trans,rest ! Restart or continue transient analysis 
*endif 
time, CURRTIME ! set current time for ANSYS solution 
deltim, TIMESTEP ! set delta-t 
£DISPNODE,fxXPFORCFUN(ICOUNT)  'apply force to DISPNODE node 
d,PV,volt,CS ! apply CS voltage to PV node 
allsel,all ! select everything to prepare for solve 
solve 
finish 
*1f ICOUNT, gt, 100,then ! turns on control] at 100th step 
G3=-11.0 
*endif 
/postl ! enter post processor 
*get,RS,node,SV, volt — get the value of the sensor voltage 
RSCheck(ICOUNT)=RS ! saves current RS in table 


* set, NODE41 Y(ICOUNT),node,DISPNODE,u,x ! get the tip displacement 
finish 


| 6 2 K k k k k k Kk kK K K K KKK K K K K K K K KKK K K K K K K K K K K K K Æ K K K K K K K K K 
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! * Control Law implementation, for next time step 
ПОДА А LUE o EE А А АД ER Seo As Е O 
*set INT,OLDINT+(RS+OLDRS/2*TIMESTEP ! approx the integral 
*set,FINT,0.95*OLDFINT-INT-OLDINT ! digital high pass filter 
*set, DBLINT,OLDDBL+H{FINT+OLDFINT)/2*TIMESTEP 
INTcheck(ICOUNT)=INT 
FINTchk(ICOUNT)=FINT 
DBLcheck(ICOUNT)=DBLINT 
OLDRS=RS 
OLDINT=INT 
OLDFINT=FINT 
OLDDBL=DBLINT 
CS=(G1*FINT+G2*DBLINT)*G3 !'CS=RS*(G1/st+g2/s*2)*G3 
CSCheck ICOUNT)=CS 
*enddo !ends the loop, ready for the next time step 


! Plot the results to screen. . 


/erase 
*vplot,, NODE41 Y(1) 


| k k k k k k k kK de ole af ole oe ole ade ol k Kk K K K K k K de a ole k K ale ale ade ol de af K k de ak ole ak K Æ K K K k k 


! * Write "NODE41Y" results to text file "truss.out" 
Y Ade k ade k k k k K k kk k k kK k ade de ol Æ K Æ K k K K k k k kk K k k k k K K K ade ade ol Æ K Æ K ak ak 


!*cfopen,piezoS_tip_c132,out, 

!*vwrite, NODE41 Y(1), 

1(E10.3) - 
I*cfclos 
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APPENDIX H. ACTIVE ELEMENT INTEGRATION DETAILS 


This appendix contains excepts from the log file used in the creation of the NPS 


space truss active model. The commands have been commented to provide clarity. All 


commands that had no relevance on the model creation were removed (for example, 


display commands) 


/PREP7 

! create element for SOLIDS 
ET,3,SOLID5 

! set parameters 
KEYOPT,3,1,0 
KEYOPT,3,2,0 

КЕҮОРТ 5,3,0 
KEYOPT,3,5,0 

! unmesh line and delete 
LCLEAR, 101 

GPLOT 

ЕЮ ЕЕЕ 0..1 

! Define Stee] material 
UIMP,4,EX, , ,200000000000, 
UIMP,4,DENS, , ‚7920, 
UIMP,4,NUXY, , ,.3, 
UIMP,4,MURX, , ,1, 
UIMP,4,PERX, , ,1, 

! place the WP at node 8 
KWPAVE, 8 

! rotate to strut axis 
wpro,-45.000000,, 
wpro,,,45.000000 
wpro,,,45.000000 

! create keypoints for the elements 
K, ,0.037545,,, 

K, ,0.189145,,, 

К, ‚0.199145 ,,, 

К, 0.2324... 

К. 

KWPAVE, 2 

К, ,-.037545,,, 

! create lines to connect keypoints 
СЭ: 8, 758 

LSIR 53, 54 

БӘЛЕ. 55 709 

ESIR. 227 0% 

LSTR, 58, 2 
KWPAVE, 54 

! shift WP for 3-axis alignment 
wpro,,,45.000000 


wpro,,,45.000000 

! create the local coord system 
CSWPLA,11,10 

! create the volume for the sensor 
BLC4,-.005,-.005,.01,.01,.01 
WESIYER 0 

KWPAVE, 56 

! create the volume for piezo 
BLC4,-0.005,-.005,.01,.01,.043 
! defines real constant steel beams 
R,8,.496e-4,.195e-9,.195e-9,.004,.004, , 
RMORES АР 

! select new lines mesh as 1 element 
! alumiumum lines 

FLST.S2 4,ORDE 2 
FITEM,5,101 

FITEM, 5,165 

CM, Y,LINE 

ЕВ. Р51 

n 

СМ. ҮЧЕ 

CMSEL,, Y 

ЕВЕ ЗВ I 
CMDEL, Y 

CMD ER YI 

! steel lines 

FLST,5,2,4, ORDE,2 
FITEM,5,162 

FITEM, 5,164 

CM, Y LINE 

БЕ Р 

pe 

CM. YLLINE 

CMSEL Y 

LESEZE S v0.2 IL. 
CMDEL, Y · 

CMDEL, Ү1 

! connector piece 
FLST,5,1,4,ORDE, 1 
FITEM,5,163 

CM, Y,LINE 


ESE „РО 

|ж 

CM, Y1,LINE 
СМЕЕ, Үү 
БЕЛЕ У 1,1, 
CMDEL, Y 
CMDEL Y1 

(m 

FLST,5,24,4, 0RDE,2 
FITEM,5,166 
FITEM,5,-189 

CM, Y,LINE 
ЕР 

|ж 

СМЕ УМЕ 
СМІБЕБ dY 
КЕБЕЖЕ 
CMDEL Y 
СЮЕ 

! mesh new lines as appropriate 
' solid alumiunum 
ТҮРЕ, 1 

МАТ, 1 

REAL, 8 

ESYS, 0 
SECNUM, , 

' mesh lines 
FEST. 2:54 ORDE.53 
FITEM.2,101 
FITEM,2,162 
FITEM,2,-165 
LMESE P5SIX 

! mesh the sensor and piezo 
ТҮРЕ, 3 

МАТ, 2 

REAL, 1 

т | 
SECNUM, , 

' mesh as one hexahedron 
MSHAPE,0,3D 
MSHKEY, 1 

ЈЕ 

FLST,5,2,6, ORDE,2 
ЕІТЕМ,5,1 
ЕТТЕМ,5,-2 

CM, Y,VOLU 
МӘБ ЕЗІ 
CM, Y1,VOLU 
CHKMSH,'VOLU' 
CMSEL STY 

|ж 

ММЕЗН, Ү1 

|ж 


CMDEL, Y 


130 


CMDEL, Y1 

CMDEIDWY 2 

! define the rigid regions at ends 
CERIGESISTALT EZ 
СЕНСОР АЕ 
CERIGPSIX ALT 
CHRIG,PSIX, ‚ALT, ss 

! change piezo to PZT material 
EMODIF,P51X,MAT,, 

! define material for quartz 
UIMP,3,EX,EY,EZ,100000000000, 
100000000000, 100000000000, 
UIMP,3,DENS, , ,23000, 
UIMP. KXX KYY KZZ IEL 
UIMP,3,MURX,MURY,MURZ,1,1,1, 
UIMP,3,RSVX,RSVY,RSVZ,1,1,1, 
UIMP,3,PERX,PERY ,PERZ 452452253 
! define material for PZT 
СІМР,3,ЕХ,ЕҮ,Е7,64.5е9, 64.5е9, 64.5е9, 
UIMP,3,DENS, , ,7600, 
UIMP,3,KXX,KYY,KZZ,1,1,1, 
UIMP,3,MURX,MURY,MURZ,1,1,1, 
UIMP;3,RSVXRSVYRSVZI 2 
UIMP,3,PERX,PERY,PERZ,1730,1730,1700, 
! activate PIEZ for PZT 

ТВЕРИ 

TBMODIF,3,3,7100 

! activate PIEZ for quartz 

ВЕ 7 

TBMODIF 3 ,3,8756 

! finished getting the control elements done 
EMODIF,P51X,REAL,8, 

! fix the elements on the active strut to material 4 
FLST 23.2 ORDE? 

FITEM,2,671 

FITEM,2,-675 
EMODIF,P51X,MAT,4, 

' unmesh strut with LPACT 
LCLEAR TTI 

! divide into two lines 

БЕЮТУ I2 T20 

FLST,5,2,4,ORDE,2 

EINEN 3127 

FITEM,5,190 

CMY. LINE 

LSE. , Paes 

'* define one element per line 
CM SY LINE 

CMSEL У 

ГЕЗ ЕСЕК: 

CMDEL, Y 

CMBET YA 

| mesh 

ТУРЕ, 1 


MAT, 1 

REAL, 1 

ESYS, 0 
SECNUM, , 

Г% 
FLST,2,2,4,ORDE,2 
FITEM,2,127 
FITEM,2,190 
LMESH,P51X 

! define mass of LPACT 
R,9,2.224, 

! mesh keypoint 
ТҮРЕ ИО 

МАТ, 1 

REAL, 9 

ESYS, 0 
SECNUM, , 


po 
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KMESH, 75 

! couple voltage DOF's at SOLIDS's 
FLST,4,4,1,ORDE,2 
FITEM,4,444 
ЕТТЕМ,4,-447 

CP. УВ Ears 
FLST,4,4,1,ORDE,2 
FITEM,4,452 
FITEM,4,-455 

CP2 VOLIPSIX 

FINISH 

/SOLU 

! define damping cónstants 
ALPHAD,9.1, 

BETAD, 1.1e-3, 
DMPRAT,0, 


! modified truss complete 
| * 





APPENDIX I. ACT_TRUSS.INP 


| ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЕЖЖ Ж 


! * act truss 3.inp x 


! * used in getting the fresh truss on line 
! * Truss Control APDL Program for ANSYS version 5.5 
! * For simple sinusoidal disturbances г 
! * Applied to the NPS Space Truss Active Controlled Model * 
! * Written by LT Carey M. Pantling 1 
! * Last Modified 14 Oct 1999 ы 


| ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖҰЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ ЖЖЖ 


ж 


! First Load Truss model with mesh and BC's 


T KKK K K KK KK K K K K K K K K K K K 


! * Define Variables 6 
| ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ 


eplot | ! gives something to look at while waiting 


! *set command will define and clear variables 


ser EREOSUETS ! disturbance frequency 

*set,CS,0.0 ! control signal to piezo, initial zero 
*set,PV,68 ! Piezo voltage node# for application 
*set,RS, ! read signal from piezo, defined only 
*set,SV,60 ! Sensor voltage node# for detection 
sen 300 ! Control gain constants, see below 
*set,G2,100 

*set,G3,0.0 

*set,OLDRS,0.0 - ! the old RS, IC 0.0 


*set, OLDINT,0.0 ! old integral, IC 0.0 
*set,OLDFINT,0.0  !old filtered integral 

* set OLDDBL, 0.0 ! prev DI, IC 0.0 

*set OLDFDBL,0.0  ! old filtered DI 

*set, STRTTIME,0.0 ! Start time = 0.0 seconds 


! set to more reasonable later 
*set,FINITIME,4.0  ! Finish Time in seconds 


TIMESTEP=1/FREQ/20 ! TIMESTEP 
*set,Pl,acos(-1) 


*set DISENODE 25 ! node 26 real truss for reading the output 
*set, APLNODE,53 ! the node where the LPACT is located 
finish 


| ЖЖЖЖж ЖЖ ЖЖЖЖЖЖЖЖ ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ ЖЖЖ ЖЖЖЖЖ ЖЖЖЖЖЖЖЖЖЖЖ ЖЖЖ ЖЖЖ 


! * Set up load disturbance, at least 20 steps per period * 
e E E E E 


*set, NUMSTEP,nint(FINITIME/TIMESTEP) 


*set, FORCFUN, 

sel EIER. 

*dim,FORCFUN,array, NUMSTEP ! array for forcing function 
*dim,FFP,table, NUMSTEP ! table for plotting 


! arrays for input (indexed), tables for output(non indexed) 


! create a magnitude based upon the frequency, from the LPACT chart 
*set, MAG,0.05437* FREQ**2-1.8874*FREQ+18.2439 ! (N) force from LPACT 


! Make a table for the force disturbance 

*do,ICOUNT,1,NUMSTEP 
FORCFUN(ICOUNT)ssin(2*PI*FREQ*((ICOUNT-1)*TIMESTEP)*MAG 
FFP(ICOUNT)=sin(2*PI*FREQ*((ICOUNT-1)*TIMESTEP))* MAG 

*enddo | | 


*set, NODE26Y 
*dim,NODE26Y,tableNUMSTEP  ! creates an array to get displacement 
ser CSC heck ! tracks the control output signal 
*dim,CSCheck,table, NUMSTEP ! gets the control signal 

* set, RSCheck ! tracks the read signal 


*dim,RSCheck,table, NUMSTEP ! copies the sensed signal 


! these may be deleted, and not be required to be copied, included for error checking 
! *set,.INTcheck 

! *dim,INTcheck,table, NUMSTEP 

! *set, DBLcheck 

! *dim,DBLcheck,table,NUMSTEP 

! *set, FINTchk 

! *dim,FINTchk,table,NUMSTEP 

! *set, FDBLchk 

' *dım,FDBLchk,table, NUMSTEP 


у жкжжжжжжжжжжтжжжжжжжжжжжжжжжжжжжжжжжжжжжжж 


1 * Do loop for loading at each time step E 
| ЖЖЖЖ FFIR TFT RITTER II ET rer I БОБО 


*do,ICOUNT.1,NUMSTEP 
! Set current time (actually, end of current step) 
*set. CURRTIME,STRTTIME+ICOUNT*TIMESTEP 
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/solu ! enter the solution processor 
*if ICOUNT,eq,1,then 


antype,trans,new ! Start new transient analysis 
else 
antype,trans,rest ! Restart or continue transient analysis 
*endif 
time, CURRTIME ! set current time for ANSYS solution 
deltim, TIMESTEP ' set delta-t 
f, APLNODE, fy, FORCFUN(ICOUNT)*0.707 ! apply force to APLNODE node 
(LPACT) 
f, APLNODE, fz,F ORCFUN(UICOUNT)*0.707 ! apply force to APLNODE node 
(LPACT) 
d,PV,volt,CS ! apply CS voltage to PV node 
allsel,all ! select everything to prepare for solve 
solve 
finish 
*1fICOUNT, gt,200,then : ! turns on control at specified step 
G3=20.0 
*endif 
/postl ! enter post processor 
*oet,RS node, SV,volt | ! get sensor voltage 
RSCheck(ICOUNT)=RS ! saves current RS in table 


*get, NODE26Y(ICOUNT),node,DISPNODE,u,y ! get the tip displacement 
finish 


| ХЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЕ Ж 


! * Control Law implementation, for next time step * 
Плата а ааа ааа кк ee ee oe 


*set,INT,OLDINT+{RS+OLDRS)/2*TIMESTEP ! approx the integral 
*set, FINT,0.95*OLDFINT+INT-OLDINT ! digital high pass filter 
*set, DBLINT,OLDDBL+(FINT+OLDFINT)/2* TIMESTEP ! second integral 
*set, FDBL,0.95* OLDFDBL+DBLINT-OLDDBL ! digital high pass filter! 
! INTcheck(ICOUNT)=INT ! copies the terms for examination 
! FINTchk(ICOUNT)=FINT 
! DBLcheck(ICOUNT)-DBLINT. 
! FDBLchk(ICOUNT)-FDBL | 
OLDRS=RS | ! cycle the terms to old values 
OLDINT=INT 
OLDFINT=FINT 
OLDDBL=DBLINT 
OLDFDBL=FDBL 
CS=(G1*FINT+G2*FDBL)*G3 ! CS=RS*(G1/s+g2/s*2)*G3 


BS 


*1f,CS,gt,60,then 
CS=60 
*endif 
*if,CS,It,-60,then 
CS=-60 
*endif 
CSCheck(ICOUNT)=CS 


*enddo 


| ak ad ale ode ode ale ade ale ale ale ade ode ale ode ale ale ade ole ade ale ale ak 


! * Plot the results to screen * 
| kk ok ok ok ok ok ok ode ЖЖ Ж Ж ЖЖЖ ЖЕЖ 


/erase 
*vplot,, NODE26Y(1) 


! voltage limiter 


! save for output 


! ends the loop, ready for next time step 


| aka ok ak de ok al ole ale ale ode ale ale k ale ale ade ale ale ade K K K k ale de ade ak ale ale ale ole K ale ale K ade ade ale ade ade ade K Æ ade 


| * Write RSCheck results to text file "ATR_e#_r#.out" * 


!* and NODE26Y results to ATN. 
!* and CSCheck results to ATS. 


ж 
ж 


| Akal ak ak ok ade ode ale ade ade ole ade ade ale ade ade ade ale ade ade ade ode ole ade ade ade ade ade ode de ale ade ale ale ade de ale ade ke ok ak ade ade ok ak 


*cfopen,AR el rl,out, 
*vwrite, CS Check(1), 
(Е) 

*cfclos 
*cfopen,AN el _rl,out, 
vwrite, NODE26Y(1), 
(Е10.3) 

*cfclos 
*cfopen,AS_el_rl,out, 
vwrite,NODE26Y(1), 
(Е10.3) 

*cfclos 
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APPENDIA J. NPS_PREP.M 


%%%%%%%%%%%%%%%%%5%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
nps_prep.m 
simulink prep program by LT Carey M. Pantling\ 
places one active member in the structure 
generates the inverse matrices and so on 
variables needed: 

Apply at поде 3 аапз 

InvMass Im 

InvMass*C ime 

InvMass*K Imk 

Application MatrixAM 

Output Select OS 


oo o9 o9 o9 o9 o9 


до до до до o9 o9 o9 


run Final_root_nps, then have the following variables in workspace 


do o9 


kk global stiffness matrix 

mm global mass matrix 

nd nodal connections for elements (nel by 2) 

coord coordinate values for the number of nodes (nel by 3) 
nnode number of nodes in system 

ndof number of degrees of freedom per node 

sdof system degrees of freedom 

bcval,bcdof boundary conditions 


до 


do oo 


до до до o9 Хо o9 


%5%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


o9 


% get the application matrix for the active element 


cel=101; % element that has active strut 
% compute elemetal dimensional data 

= свека па се! 2).1)-соотачпа се Tee 
ау=ввова паее 2 ) 2) -свекасла(еве 20 
adz=coord (na (cel 2) 3)-совка(па(се1! 1) =): 
Тепда= (ах 2+ду 2+82 2)70.5; 


compute element rotation matrix 
local aligned with x along axis, y,z orthogonal 
Since tube segments, y,z do not matter 
DCM defined as {x}g=g[c]1*{x}1; g[c]1 is the transformation matrix 
ul=[(dx/leng) , (dy/leng) , (dz/leng) ] ; 
er eress MM 
mags (cr(1)^24cr(2)^24cr(3)^2)^0.5; 
if mag--0 % if {ul}x{i}-=0 then 
OA=[0,1,0;-1,0,0;0,0,1]*ul'; $ rotate z-90degrees 
else 
OA-[0,0,-1;0,1,0;1,0,0]*ul'; $ rotate y-90 degrees 


до до до 


до 


епа 
vlSeross (ul Ол) о. 
wl=cross (ul, vl) > wl=w1/ WI Zr Zul ee: 
рем= [чі Уту 5 оо е 

uld2) 2) уа 
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ие la, wi Be 
r=zeros (12); 
ze 3 33) DCM:. 
r (476, 4:6) =DCM:; 
77917 :9) ЗБЕ 
т (ШОО УО: т2)=Ю М1 
% ТЕ ЕТ 1сСд=г' 
am=zeros (12,1); 
am(1)=-1; am(7)=1;% axial transfer matrix 
am=r*am; % app matrix=gcl*atm 


ео (12). 
isom(1)=-1; isom(7)=1; 
lsomsisom*r'; 


AM=zeros (sdof,1); % Initialize the Application matrix 
ISOM=zeros (1,sdof) ; % Initialize the isolator matrix 
start= (nd (cel, 2) -1)*6+1; 

AM(start:start+5)=am(1:6); 

ISOM (start :start+5)=isom (1:6); 
г:акЕ=(га(сеј 1) DE +T;; 

AM(staret:.starcs5) =am 7: 12) 

ISOM (start :start+5)=isom(7:12); 


до 


[m] (a)« [c] (v) « Ik] (x) s (£) 

(a) «inv [m] [c] (v) *inv [m] [k] (x] 2inv[m](£) 

setting BC's sets imk values to 1.0 the imf values to BCval 
(572 [т] +5 [2] + [К] ) (х) = (2) 


до до до 


% Inverse Mass 
$ requires squishing the mass matrix to the displacement terms 
inverting, expanding 
msq=mm; $ preserves the original mass matrix 
for index=nnode:-1:1 
% remove the back terms first so as to not let the matrix be confused 
start=(index-1)*6+4; 
fini=start+2; 
ШОО Бак ты: = |. 
usqi-start:ftrini)s 
end | 
invmsq=inv (msq) ; 
% repopulate the Im matrix 
Im=zeros (sdof) ; 
for index=1:nnode 
for j=l-nnode 


start= (index-1) *6+1; ' startsq= (index-1) *3+1; 

fini=start+2; ‚ Finisq=startsq+2; 

Im(start:fini,start:fini)-invmsq(startsq:finisqd,startsg fine ge 
end 


end 


% Inverse Mass * C (no natural damping) 
Imc=zeros (sdof) ; 
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% Inverse Mass * K 
Imk=Im*kk; 


apply the boundary conditions, doesn't matter since all zero 
ШП, ЕЕ Д = еар усо (етт Е СЯО Бета). 


оо до 


remove the constrained points from the mass and stiff matrices at 
nodes 1,2,27,28 we're not tracking boundary forces, so we don't care 
% what the forces are and can remove from model 


o? oe 


sdof=sdof-24; 

% reduce mm 

Im(157:168,:)=[]; % nodes 27,28 
Im(:,157:168)=[]; 

ти (1:12 затар; $ nodes 1,2 
Ent: 1:12) 21]. 

$ reduce Imc 

Imc (157:168,:)=[]; % nodes 27,28 
ime (oon 68) = [] > 
Ітс(1:12,:)- (1; 
met 1:12) -(1- 

$ reduce kk 

Imk (157:168,:)=[]; 
Imk(:,157:168)=[]; 
зо RS A % nodes 1,2 
Imk(:,1:12)=[]; 

% reduce AM 

AM(157:168)=[]; % nodes 27,28 
AM(1:12)=[]; % nodes 1,2 

$ reduce ISOM 

ISOM(157:168)-[]; % nodes 27,28 
1Ї$ОМ(1:12)=[]; $ nodes 1,2 


оо 


nodes 1,2 


nodes 27,28 


oo 


$ Apply at node 3, (now node 1) 
aan3=zeros (sdof,1); 
aan3(1:3)=1; 


$ output select gets y translation of node 3, (now node 1) 


OS=zeros(1,sdof); 
OS(2)=1; 


12% 











APPENDIX K. HP 35665A DYNAMIC SIGNAL ANALYZER SETUP 


ң 


The HP 35665 A DSA has several different keys that are available to the user to establish 
a configuration. These keys are of two types hard keys, which do not change value and 
soft keys, which do change value as various hard keys are pressed. Hard keys will be 
represented by [Hard Key]; soft keys will be represented by [Soft Key]. АП the specifics 
on the commands used are contained in Reference 29. Hewlett-Packard, HP35665A 


Dynamic Signal Analyzer Operator 's Reference. 


Step 1: Turn on the spectrum-analyzer. 
The 1/0 (power) switch 1s on the front panel of the machine. 
Step 2: Define the analyzer instrument mode. | 
[Inst Mode] [2 CHANNEL]. 
Step 3: Define the analyzer’s frequency bandwidth. 
[Freq] | SPAN ] enter 200 Hz, 
[RESOLUTION Jenter 800 Lines, and 
[RECORD LENGTH] enter 4 Seconds. 
Step 4: Select the Display Format. 
[Disp Format] [ UPPER/LOWER |. 
Step 5: Display frequency magnitude response on trace A. 
[Meas Data] [ FREQUENCY RESPONSE |. 
[Trace Coord] [ MAGNITUDE |. 
Step 6: To Display frequency phase response on trace B. 
[Active Trace] will toggle the active trace to trace B 
[ FREQUENCY RESPONSE ]. 
Step 7: Define Trace B coordinates to phase. 
[Trace Coord] [ PHASE ]. 


Step 8: Establish Random noise Source. 
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[Source] [ LEVEL ] enter in 1.00 [ Vpk ]. 
[ RANDOM NOISE ] 
[ SOURCE ON/OFF |] toggles on the active source 
Step 9: Create Frequency Resopnse Window. 
[Window] [ HANNING |. 
Step 10: Establish Averaging for Random Signal. 
[AVG] [ VECTOR ] preserves the phase information . 
[ NUMBER ] select 20 average cycles 
Step 11: Save data on 3-1/2 inch floppy disk. 
[Disk/Util] [ SAVE DATA] 
[ SAVE TRACE ] 
[ INTO FILE ] type in filename at attached keyboard. 
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APPENDIX L. DSA_PLOT.M 


ххххххххххххххххххххххххххххххххххххххххххххххххххххх 


ao до 


ж ж 
* HP 35665A Data processing utility * 
* Designed for the NPS Space Truss Experiment 5 
* C:\Space_truss\truss99\Experiments\hp_dsa * 
x ж 
ж ж 
ж ж 
ж ж 


до о оо 


By LT Carey M. Pantling, USN 
Last Modified 8/6/99 


A A a 


ххххххххххххххкхххкхххххххххххххххххххххххххххххххх 


oe 


cd c:\space truss\truss99\experiments\hp dsa 


1осоаа а99_051а.та© t load the data file from the HP-35665A 
creates the following variables 

© 1 the vector of nlines+1 of complex results 
o2i1x0 

оз хр. frequency resolution 

Estart=0: 

ога $ print the freq resolution to WS 

fend=50; 


oe 


A оо до 


fEreq=\he tart. O2 1 lect fendis 
figure (1) 


Subplot(2 mu 3 
semilogy (freg,abs (02il)); % plots the magnitude 
axis([fstart,fend,1e-4,1]); 
xlabel('Frequency (Hz)'); 
ylabel('Magnitude'); 
title('Magnitude and Phase plots (0-50 Hz)'); 
grid; 


Gubplob(2 14217 
ph=atan2 (imag (02il),real(o0o2il)); % creates phase +/- pi 
plot (freg ph leo, pt). % plots phase in degrees 
axis([fstart, fend, -180,180]); : 
xlabel('Frequency (Hz)'); 
ylabel ('Phase (deg)'); 
grid; 
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APPENDIX M. ACQ DATA.M 


хх ХХХ ХХХ ko AX ХХ 


оо оо 


* * 
* Data Acquisition subroutine for dSPACE ж 
* Designed for the NPS Space Truss Experiment * 
* By LT Carey M. Pantling, USN * 
* Located C:\Space truss\Experiments\Experiment * 
x ж 
ж ж 
x * 


A до 90 909 


Last Modified 8/23/99 


% 


o9 oo 


ххх хх ХХХ Xx 


oo 


% verify that the correct directory is being used 
cd c:\space_truss\truss99\experiments\experiment_2 


$ select 051103 board for use with MLIB and MTRACE 
mlib('SelectBoard', 'DS1103'); 
mtrc1103('SelectBoard','DS1103'); 


$ check if the application ST controller 3.ppc is running on dSPACE DSP 
DemoApplName = lower(['c:\space truss\truss99\st_controller 3.ppc']); 
if mlib('IsApplRunning 
ApplInfo = mlib('GetApplintos 
if strcmp (DemoApp1lName, lower (Applinfo.name)) ~= 1 
err msg - sprintt( еше СЕ” пе РРС У. 
'applicationin*** ''$s'' running!', DemoApplName) ; 
сво Еве С. 
end; 
else 
err msg - sprintf('*** This m-file needs the PPC ‘,. 
‘application\n*** ''%$s',' running!',DemoApplName); 
error(err msg); 
end; 


$ select variables to be traced and obtain their descriptors 

var names =('Model Root/B:input->ADC'; 
'Model Root/B:Sum'; 
'Model Root/S:Accelerometer 1( 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 
'Model Root/S:Accelerometer 


HS 


i 


Ph и бо бол бо КУ КУ КУ 
Wh) Bt) Pp Ооо o pp 


ut me . 
e 
. . 


var = mlib('GetTrcVar',var_names); 


% set data acquisition to be performed by service number 1 (default) 
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=т= о; % sample time (sec) 
ES-500, % sample Frequency (Hz), set by simulink parameters 
Пете. t€ number of samples 


mtrcl103('Set','TraceVars', var, 'NumSamples',NS); 


mErcIIOS3Q'Startcapture е % stare capeune оп 05 1703 
while mtrc1103 ('CaptureState')-=0,end % wait until capture is done 


% Fetch after capture is complete 
out data = mtrc1103('FetchData') ; 


% Plot code extensively modified from 'Graph.m' 
% originally written by LT John Vlattas and LT Scott Johnson 5/10/98 


timevec=[1/FS:ST/NS:ST]; 
fOr = 4 


eval (['Y' num2str (n) n =Soüut datali, :); 1); 
% Separates Data to vectors 
їп ез 
figure (1) 
subplot (2,1,n) 
rM 


plot (timevec,eval(['Y' num2str(n)])) 
title('PCB Force Sensor - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mv) ') 
elseif n == 2 
plot (timevec,eval(['Y' num2str(n)])) 
title('Output of Controller - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mV) ') 
end 
orient tall 
ешселте те 6 
figure (2) 
subplome 1 n2» 
ilt us | 
plot (timevec, eval Y saruncstrin) ) m) 
title('Node 41 - X-Axis - Time Data') 
xlabel('Time (msec) ') 
ylabel ('Amplitude (mv) ') 
elseif n == 
plot trimevec,evale@l у пио 55 (п) 1)" ) 
title('Node 41 - Y-Axis - Time Data!) 
xlabel ('Time (msec)') 
vlabel ('Amplitude (mV) ') 
elseit — 
plot (timevec,eval(["Y” Mu2str(n))),'m") 
title('Node 41 - Z-Axis - Time Data') 
xlabel('Time (msec)') 
ylabel ('Amplitude (mV) ') 
end 
orient tall 
elseif n < 9 
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figure (3) 
subplot (3, 1,sGieS)) 
if n == 
plot (timevec, eval. ([{'Y' @num2stmin) ]), ш") 
title('Node 18 - X-Axis - Time Data') 
xlabel ('Time (msec) ') 
ylabel ('Amplitude (mV) ') 
elseif n == 
plot (timevec,eval(['Y' num2str (n)]),'m') 
title('Node 18 - Y-Axis - Time Data') 
xlabel ('Time (msec)') 
ylabel ('Amplitude (mV) ') 
elseif n == 8 
plot (timevec,eval(['Y' num2str(n)]),'m') 
title('Node 18 - Z-Axis - Time Data') 
xlabel('Time (msec)') 
ylabel ('Amplitude (mV) ') 
end | 
orient cari 
elseif n < 12 
figure (4) 
subalsel(>37171.(n-8)) 
if n == 
plot (timevec,eval(['Y' num2str(n)]),'m') 
title('Node 49 - X-Axis - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mV) ') 
elseif n == 10 
plot (timevee, eval (u num2str(n)]); 'm' ) 
title('Node 49 - Y-Axis - Time Data') 
xlabel ('Time (msec)') 
ylabel('Amplitude (mV) ') 
elseif n == 11 а 
plot (timevec,eval(['Y' num2str(n)]),'m') 
title('Node 49 - Z-Axis - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mV) ') 
end 
orient tall 
elseif n < 15 


figure (5) 
subplot (3,1, (n-11)) 
if n == 12 


plot (timevec,eval(['Y' num2str(n)]),'m') 
title('Node 26 - X-Axis - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mV) ') 

elseif n == 13 
plot (timevec,eval (['Y' num2str(n)]),'m') 
title('Node 26 - Y-Axis - Time Data') 
xlabel('Time (msec) ') 
ylabel('Amplitude (mV)') 

elseif n == 14 
plot (timevec,eval(['Y'! поо Е а) 
title('Node 26 - Z-Axis - Time Data') 
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xlabel('Time (msec) ') 
ylabel('Amplitude (ту)') 
end 
orient tall 
end 
end 


% Takes Power Spectral Density of Given Vectors For Comparison 
Pxx1=fft (100% vw 

Pxxl- (real(Pxx1).^2«imag(Pxx1).^2).^0.5; 

Pxx1((NS/2+1) :NS)=[]; % clears the redundant upper portion 
F1=FS/NS* [0:NS/2-1]; % vector for plotting the fregs 


% Load the uncontrolled case 


load Pxx2 % PSD vector 

€ Finds Maximum Value of the T For Trial In Question 

u = max(10*log10 (Pxx2 (10:NS/2)) 

¢ = max (10*10g10(Pxx2 (10-NS/ 2) | 

diff = u - c $ dB reduction due to control 


define trial number 
$ trial=6;°  $ trial number located in workspace 


до 


figure (6) 

% plot (F1(1:400),10*1og 10 Ри A090) 

plot(F1(1:400),10*10g10(Pxx2(1:400)),'r.', 
F1(1:400) ,10*1log10 (Pxx1 ( een) ID 1) 

Gia 

€ title(['Power Spectral Density - Uncontrolled']); 

title(['Power Spectral Density-Controlled vs Uncontrolled (Trial ', 

пишер ЕР)" 

ylabel('Power Spectrum Magnitude — JE 

xlabel('Frequency (Hz)'); 

legend ('Uncontrolled','Controlled'); 

% save trial information in appropriate file 

eval (['save trial',num2str(trial),” Y* Pxxl FE as 
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APPENDIX N. DATA_PROC.M 


% Ххх 4544445424 44 4 4 444 4 4 4 4 d dd dd 44 4 

$ * Data Processing routine for dSPACE data * 

% * Designed for the NPS Space Truss Experiment * 

$ * By LT Carey M. Pantling, USN * 

% * Last Modified 10/16/99 * 

% ХхХххХ хх RAR AAA AAA A AAA A AAA AAA AA AAA AAA AAA RARA kxk 

ТЕ. + sample time (sec) 

E5=500; + sample Frequency, (dz), set by simulink parameters 
NS=ST*FS; < number of samples 


Plot code extensively modified from 'Graph.m' 
originally written by LT John Vlattas and LT Scott Johnson 5/10/98 


до до 


% trial number stored as variable in workspace 
% load trial information from appropriate file 
eyal( ("load tria PP Fnum2sftr (Erialy Ух Ру јан"): 


timevec= [1/FS:ST/NS:ST]; 
fOr n = 22а 
ШЫР лез 
figure (1) 
sübplot (2, mi) 
Dn 
plot (timevec,eval(['Y' num2str (n)])) 
title('PCB Force Sensor - Time Data') 
xlabel ('Time (sec)') 
ylabel('Amplitude (mV) ') 
elseif n == 
plot (timevec,eval (['Y' num2str (n)])) 
title ('Output of Controller - Time Data') 
xlabel ('Time (sec)') 
ylabel('Amplitude (mV) ') 
end 
orient tall 
elseif n< 6 
figure (2) 
subplot (3 1, (1-27) 
if n == 
Plot (timevee ела Fnum2ser(n)]) то) 
title('Node 41 - X,Y,Z-Axis - Time Data') 
ylabel ('Amplitude (mV) ') 
elseif n == 
plot (timevec,eval(['Y' num2str(n)]),'m') 
ylabel ('Amplitude (mV) ') 
elseif n == 5 
Plot tEimeveec,eval(?Y num2striny it. om) 
xlabel('Time (sec) ') 
ylabel ('Amplitude (mV) ’) 
end 
orient tall 
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elseif n < 9 
figuret3) 
subplot (27.1, (5579 
Реп. 6 
plotitimeves, eval У аи (a) D mds) 
title('Node 18 - X,Y,Z-Axis - Time Data') 
ylabel('Amplitude (mV) ') 
elseif n == 7 
plot (timevec,eval(['Y' num2str(n)]),'m') 
ylabel('Amplitude (mV) ') 
elseif n == 8 
plot (timevecheval( (Nine ste (a) 19,0% 
xlabel ('Time (sec)') 
ylabel('Amplitude (mV) ') 
end 
orient tall 
elseif n < 12 | 
figure (4) 
subplot (3,1, (ae) ) 
ifn = 
рјов је те“есема Ра ита Ела Mm.) 
· 12 Је ('Моде 49 - X,Y,Z-Axis - Time Data') 
ylabel ('Amplitude (mV) ') 
elseif n == 10 
plot (timevec ела Ти ани em 
ylabel('Amplitude (ту)') 
elseif n == 11 E 
plot (timevec,eval(['Y' num2str(n)]), 'm') 
xlabel ('Time (sec) ') 
ylabel ('Amplitude (mV) ') 
end 
orient tall 
elseif n < 15 


figure (5) 
SUDPLOE тит 
if n == 12 


plot (ermevec eva” num2stri(nu)]) m?) 
title('Node 26 - X,Y,Z-Axis - Time Data') 
ylabel('Amplitude (mV) ') 

elseif n == 13 
plot(timevec,eval(['Y' nuestra), mo) 
ylabel ('Amplitude (mV) ') 

elseif n == 14 

plot (timevec,eval(['Y' num2str(n)]),'m') 

xlabel ('Time (sec)') 

ylabel ('Amplitude (mV) ') 


end 
Оттеп Ба 
end 
end 
[Pxx1,F1]=PSD(Y13,2048,FS); % creates PSD of node 26-Y 


о 


% Load the uncontrolled case 
load Pxx2 Pxx2 
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6 Finds Maximum Value of the Power Spectral Density For Trial In 
Question 

п -йтах (20510010 (Рхх2 (10:100))) 

с тах (20*10910 (Pxx1l(l0:100))) 

атт y= MC $ dB reduction due to control 


figure (6) 

% plot (F1(1:100) ,20*log10 (Pxx1(1:100))); 

Plot (FRE ОВОО ВИ P 22910900, !r 77,27 (12100), 
202410910 (Pxx1(1:100)), D=) 

grid; 

% title(['Power Spectral Density - Uncontrolled']); 

title(['Power Spectral Density - Controlled vs Uncontrolled (Trial ', 

: num2stritr al). p 

ylabel('Power Spectrum Magnitude (dB)'); 

xlabel('Frequency (Hz) '); 

ss legend (* Uncontrolled 0); 

legend *Uncontrolled', "Controlled" 0); 
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APPENDIX O. EXPERIMENTAL RESULTS 


- 


EXPERIMENT 1 12.50 Hz 1.2 Vpp 


0 | 0 | 0 | 0 | 46:728 | - | Uncontrolled у 
ENERO — ЛЛ -- "tcl | | 
_2 | 350 | 0 | 2.00 | -57.6053 | 10.8926 | Long-termunstable | 
3 | 350-і шайы. A S | 
Ca ¡OI | 
| 6 | 300 | 2500 | 1.75 | -62.8023 | 16.0895 |  UnstableinlOsec | - | 
|7 |300 | 1000 | 1.75 | -61.0617 | 143489 | Unstableinl5sec_| —— 
| 8 | 300 | 500 | 1.75 | -62.6263 | 15.9135 | Unstablein20sec_ | 
_9 | 300 | 200 | 175 | -60.7286 | 14.0158 | Borderlineunstable | Y | 
| 10 | 300 | 0 | 1.75 | -623806 | 15.6678 | Changed AirTank | | 
ПО ССОО ЕЕ О И Ис И NE 
[12 | 300 | 100 | 2.00 | -62.5081 | 15.7953 | Borderline stable | 
|.13 | 300 | -100 | 1.75 | -64.3968 | 17.6840 | 04/25 
| 14 | 300 | -200 | 1.75 | Lost | Los | DataFilelost | | 
| 15 | 300 | -400 | L.75 | -65.2538 | 18.5410 | BestCase, this Trial | | 
| 16 | 3000500 | 1.75 ANA A ы 8886 
| 17 | 300 O A rara 07290. ыы | __ 
_18 | 0 | 400 1175| Lost | Lost |  DataFileLost | | 
ШЕШЕ! (77777: | — 
20 | CATIA AIM | 
AA | ОЧ ИТ ШЕКЕЛ ж A o ШЕ 
_ 22 | 0 | -4000 | 1.75 | -49.3786 | 2.6658 | Unstabledsec | | 
_ 23 | 300 Ето шта сес er A 
24 | 0 | 0 | O0 |-72137 |  - __| Uncontrolled 1.0Vpp | — 
| 26 | 300 | -200 | 2.25 | -74.9435 | 2.8058 | Borderline stable | | 
| 27 | 300 | -200 | 2.50 | -73.5656 | 1.4279 |  Borderinestable | 
_ 28 | 300 | -200 | 3.00 | -70.9681 | -1.1696 |- Unstable | | 


Table 19. Experiment 1 Results 
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Figure 45. Exp.1 Trial 0 Controller and Node 41 Response 
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Figure 46. Exp. 1 Trial 0 Node 18 and 49 Response 
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Figure 47. Exp. 1 Trial 0 Node 26 Response and PSD 
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TRIAL 5: IGain=300, IGain=5000, SG=1.75 


PCB Force Sensor - Time Data 


| 
ТО, hy T h 


0.04 | 
0.03 | | | 


l 


Amplitude (mV) 
о 
( c. 
N 





0 2 4 6 8 10 12 14 ^6 18 20 
Time (sec) 


Output of Controller - Time Data 


Amplitude (mV) 









0 2 4 6 8 10 12 14 16 18 20 
Time (sec) 


Node 41 - X,Y Z-Axis - Time Data 


| 
„ыш 
A 


о 
о 
— 


< 0.005 E | 
Е |! j 
© а Sl bans думан ша ea tilt AI guy Kutgil Kay » aa! д, РА е7) | ые à Т Тр! 
5 ОМ МИ И млевено б: i 4 {е ~ 4 шш айнам Wit thm Ae ү 
PC | 
& .0.005 

-0.01 

0 2 4 6 8 10 12 14 16 18 20 


| 
| 
{ pie vil |: Y at + ie | | hg | 


M (ип Ити а. fi 
| ap lu Asie ne bt к PUR FCR PRI і | 


mi ШҮН „neh 


; Mi ith " ih » te ty! iM 


Amplitude (mV) 





0 01 












S 0.005 | 

Е ШИШ Г 

= Иш M ШҮ D VE "n Мың y % vindi TT fin үң ET |, si 
Ф Nu Ing) vi У Y ч) ТИ | | 

E d МІ itil ое ЙЫ. ШТ НЙ T ү " ah d Ж | ү ч Ш N Hi int 5 yea P 
с $, 

Е i ' } 


0 2 4 6 8 10 12 14 16 18 20 
Tıme (sec) 


Figure 48. Exp. 1 Trial 5 Controller and Node 41 Response 
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Figure 49. Exp. 1 Trial 5 Node 18 and 49 Response 
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Figure 50. Exp. 1 Tnal 5 Node 26 Response and PSD 
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TRIAL 9: IGain=300, IIGain=200, SG=1.75 
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Figure 51. Exp. 1 Tnal 9 Controller and Node 41 Response 
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Node 18 - X. Y,Z-Axis - Time Data 
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Figure 52. Exp. 1 Tnal 9 Node 18 and 49 Response 
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Node 26 - X,Y Z-Axis - Time Data 
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Figure 53. Exp. 1 Trial 9 Node 26 Response and PSD 
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EXPERIMENTE 2> 13.81 Hz] 2Vpp 






| 0| 0 ' O | о | 505009 | - | Uncontroled | 
1 |30 |0 | 200|..65.1144 | 1453356 ШЫ | 
| 2 | | 0 |200! -67.6224 | 17.0415 | Unstable | 1 
[3 | 350 | O 11.75 | -68.8184 | 18.2484 | Borderline unstable | — 
_4 | 300 | 100 | 1.75 | -68.8293 | 18.2375 | Borderline stable | | 
__5 | 300 | 200 | 1.75 | -70.4061 | 19.8251 | Borderofstability | | 
io | ОШ ЫШ пошто ж” | 
27 |300 | -200 | 1.75 | -73.3524 | 22.7714 | NotAbletoDuplicate | — 
__8 | 300 | -300 | 1.75 | -69.6011 | 19.0202 | BESTOVERALL | У 
[9 | 300 | -200 | 2.00 | -65.7410 | 15.1601 | CheckofTia7 | — 
[10 |30 | 0 | 2.00 | -66.6074 | 16.0264 | (| - 
_ ипи | 0 | 200 | 2.00] -500963 = 
EAE 2 a. | — 
сз | ПИШ ВЕ ЕЕ | __ 
[14 | 0 | 5000 | 2.00 | -50.5151 | -.0658 |  Unstable9see__ | 
SAO ООО БО ЕТС ЕС Өө В 
ООН ty OC ОЕ И 
[17 | 0 | -5000 | 2.00 | -51.5259 | 0.9450 | UnstableSsec_ | 
[18 | 0 | O | 0 | -84325 | - | Uncontrolled10Vpp | — 
| 19 | |i cee | 
_ 20 | 300 | 0 1225| -98.5186 | 14.1761 | NoisySignal | 1 
| 21 | 300 | 0 | 3.5 | -96.4280 | 12.0855 | _Unstable2sec | | 


Table 20. Experiment 2 Results 
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Figure 54. Exp. 2 Trial 8 Controller and Node 41 Response 
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Figure 55. Exp. 2 Trial 8 Node 18 and 49 Response 
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Node 26 - X, Y.Z-Axis - Time Data 
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Figure 56. Exp. 2 Tnal 9 Node 26 Response and PSD 
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EXPERIMENT 3 16.75 Hz 1.2 Vpp 


HOLEN | .0 5280€ 2 MEL Uncongoled — |" — 
_1 | 300 | 0 | 2.00 | -68.7235 | 15.9385 | Borderline Unstable | | 
_2 | 350 |707 27200 |: 70908109] 9018/35422 0 A stable | | 
[3 | 350 | 0 11.75 | -672985 | 145135 |  SlowlyUnstable | | 
г 4 773006 ONIEEN Oc. ee ee 10. РР 8 С 
[5 | 300 | 200 | 1.75 | -70.4232 | 17.6382 | Borderline Unstable | | 
6 | "ORO о у | | 
7 [0 | -200 | 1.75 | -57.2934 | 4.5084 | Noisy Output |: 
[8 | 300 | -100 | 1.75 | -70.7148 | 17.9299 | BestCase, This Trial | | 
- 9 | 300 | -200 | 1.75 | -70.7771 | 179921 | Borderline Stable | —— 
[ 10 | 3002/25007 753 705380 (20805070 9 004 P N 
[11 | 300 | -200 | 2.00 | -74.3528 | 21568) Unstable | 
| 12 | 0 | -200 | 2.00 | -57.2314 | 44465 | № | | 
| 13 | 0 1 ЕЕ о WI E Posa | 
—14 | © | -1007/ 2001 -56214% | "34298-19 %Ұ%% 5000 3| — 
| 15 | 0 | 0 |0 j| 565548 |  - | Checksuncontrolled_| | 
SE Er | озы мН 


Table 21. Experiment 3 Results 
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EXPERIMENT 4 12.50 and 13.81 Hz, IGain=300, HGain=100, SG=1.75 
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__41 | 180 | -73.9794 | 7.8041 | — 
ОЕ” |Ж” СӨЗБЕН Е ИН 
co о Ж. 
[SATA AA л 


Table 22. Experiment 4 Results 
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TRIAL 51: Phase = 240 degrees 
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Figure 57. Exp. 4 Trial 51 Controller and Node 41 Response 
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3 Node 18 - X,Y Z-Axs - Time Data 
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Figure 58. Exp. 4 Trial 51 Node 18 and 49 Response 
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Power Spectral Density - Controlled vs Uncontrolled (Tnat 51) 
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Figure 59. Exp. 4 Trial 52 Node 26 Response and PSD 
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Power Spectrum Magnitude (dB) 


EXPERIMENT 5 Random Noise 1.0 Vpp 10-35 Hz window 
IGain=300, IIGain=100, SG=1.75 


Be En О EEE 
BI ne 
AE 
| Mean | -77.65034/-5 | | -84.3384 4-5 | 
Table 23. Experiment 5 Results 
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Figure 60. Exp. 5 Average Power Spectral Density 


APPENDIX P. PROC_ANSYS.M 


ХхххххххххкхххХхх лс ХХХАХХХХХХХХАХХХАХХХАХХАХХХХХХХХАХХХХ 


оо 


$ * Data Processing subroutine for ANSYS * 
% * Designed for the NPS Space Truss Experiment * 
% * By LT Carey M. Pantling, USN * 
% * Last Modified 10/21/99 * 


хххххкхххх хх ХКК ХХХ ХХХ RR 


до 


freq—ii- 75. 3 driving frequency 
FS=20*freg; % sampling frequency 
NS=940; % number of samples 
time=[1/FS:1/FS:NS/FS]; $ time vector 
expt=1; % experiment number 


% АМ - node 26-y output 
% AR - output to piezo 


$ AS - Input from force transducer 
tor - 28 % runs through all the data 
% load the data and get into a good name 
ШЕ то $ gives leading zero 
eval(['load AN e',num2str(expt),' O',num2str(n),'.out;']); 
eval(['AN- AN et, num2stfr (ере), Соате (а: 
eval(["load AR’e7-num2str\expt), ' 0T nuUm2str(n) em is 
eval(['AR = AR e',num2str(expt),' O', Dumnesstrp XD | 
evalií[ load. AS: пиш? expt), 0 num2strin "ont |), 
eval (['AS = AS e',num2Str(expt),' 0T7,num2strm),, "225 
else % no leading zero 
eval (['load AN e! num2str (expt), ' M num2str(n),'.ouťt; m); 
eval(['AN = AN e", num2str(expt),'." num2stre(n),';'1); 
eval(['load AR e',numa2str(expt),' S'Snmnum2strín),'.out;']); 
eval (['AR = AR e' num2str(expt),'" w"enüum2str(n);';*']1Y; 
eval([|'wogd9AS e',nuüum2stm c pt АШ 5 ЕЈ (ај о ОНЕ lOs 
eval(l AI =7257eY, num2sernespe) alme sean). У 
епа 
Figure): % plot the raw data 


subploe >, 1.1), 
plot (time, AS); 
title('PCB Force Sensor and Controller Output') 
ylabel('Amplitude (V)'); 
ах15([0 4 -1 11); 
subplot? i, 2); 
plot (time,AR); 
утаре ү уусу; 
ахіс((0 4 -80 801); 
subplot (3,1,3); 
plot (time,AN); 
title('Node 26- -y Motion') 
ylabel (' (meters) '); 
xlabel('Time (sec)'); 
eval (['print -djpeg t, num2Str(n),'figl IP J); 
figure (2); * plot the processed data 
[Pxx1,F1]=PSD(AN,512,FS); $ creates PSD of node 26-Y 
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if n== g% plot only the uncontrolled PSD 


plot (F1 (14100) 2907 10910 (Pxx1 (1: 100))?); 
title(['Power Spectral Density - Uncontrolled']); 
ylabel('Power Spectrum Magnitude (dB)'); 
xlabel('Frequency (Hz)'); 
legend ('Uncontrolled',0); 
axis([0 40 -300 -100]); 


grid, 
PXXZ= EXI 
save Pxx2 Pxx2; % save the uncontrolled PSD 
else $ evaluate the controlled cases 
load'Pxx2; % load uncontrolled case 


plot (Fi ( Bor 20° Tograomp 22 1(12100)), ?r 779700) 
20* loglO(Pxxl (1: 10058) b--") 
title(['Power Spectral Density - Controlled vs Uncontrolled', 
' - Trial numeserm).')'])- 
ylabel('Power Spectrum Magnitude (dB) '); 
xlabel('Frequency (Hz) '); 
legend ( 'Uncontrolled','Controlled',0); 
axis([0 40 -300 -100]); 
атша ; 
eval (["print -dgpeg = саме, Ета ра I; 
с = пах(рос од ет на oe 
ve = max (20410010 (Pxx] (102100 oe 
ПЕТ “кш % dB reduction due Со соп го! 
үебпігсеіп,:)-|(п u c diff] 


end % of plot routine 
end % ОЕ loop for rrrat ante 
results % spit results to screen 
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APPENDIX Q. ANSYS RESULTS 


Series 1: 11.30 Hz 


01010 |0 | -103.5685 | - | Uncontrolled Case | _ 
MEA. e o | AA. ASE. a 
ATA E AIR жәнен O AA 
з | о ы AAA | 
[| 4 ПЕРНЕ E ë ae 
ls JEEP A ASS ЕЕЕ | 
[6 | 300 | 0 |-125 | -115.3466 | 11.7781 |  ClosetoSaturated | - 
[7 | 300 | 0 |-150 | -115.3658 | 11.7973 | _Saturation Limit | | 
[8 | 300 | 0 |-175 | -115.4091 | 11806 | Saturated | 
[9 | 300 | 0 | -200 | -115.4285 | 11.8599 | SaturatedIGain | Y | 
| 10 | о | 50k | -150 | -127.0601 | 23.4915 | _Saturated Gain | Y | 
| 11 | 0 | 25k |-150 | -121.8247 | 18.2562 | Saturated | | 
12 | 0 | 20k |-150 | -123.5433 | 19.9748 | Saturated — | — 
[13| 0 | 15 | -150 | -125.1943 | 21.6258 | Saturated | 
14 | 0 | 10К | -150 | -128.5719 | 25.0034 | Saturated | 
| 15 | 0 | Sk |-150|-129.6683 | 26.0998 | Slightly Saturated | Y 
| 16 | O0 | 5К |-150 | -126.0206 | 22.4521 | Везё ОуегаН Сазе | У_ 
О О О В О А | 
SAN SA ыл ы т шшш NNI 
о ое КИ БЕРЕ И ВОО ш 
AE AE ЕЕ AAA A аш 
IIA EUA AA И A 
DAA ШЕЛ E ОВ 
| 23 | 300 | 2k | -75 |-105.9514 | 2.3829 |  Questionable | | 
_ 24 | 300 | 4k | -75 | -110.5684 | 6.9999 | Questionable | _ 
SS TEM MS AAN | Б НА 
26 econ a A 
C27 e EAN AE И 


Table 24. ANSYS Series 1 Results 
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Tnal 9: IGain=300, IIGain=0, SG= -200. 


PCB Force Sensor and Controller Output 
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Figure 61. Senes 1 Trial 9 Results and PSD 
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Trial 10: IGain-0, IIGain-50,000, SG- -150. 
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Figure 62. Senes 
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Trial 15: IGain=0, IIGain=8,000, SG= -150. 


PCB Force Sensor and Controller Output 
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Figure 63. Series 1 Tnal 15 Results and PSD 
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Trial 16: IGain-300, IIGain-5,000, SG- -150. 


PCB Force Sensor and Controller Output 
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Figure 64. Series 1 Trial 16 Results and PSD 
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Series 2: 11.75 Hz 


0 | o | o | o |-1046728 | | Uncontrolled Case | —— 
_1_| 360 Оо ш ШО о 
ШЕШЕГІ” о МАРУ ee асс) ж” 
o o ШКО ОЕ к ИИ 
_4_| 500" у пш к ши к ИВ 
| 5 | 300 | 0 | | -119,2355 | та 62 | 71 ЕЕ 
[6 | 300 | 0 | -125 | -119.5663 | 14.8936 | — Saturatedlimit — | — 
| 7 | 300 | 0 | -150 | -119.6658 | 14.9930 | Saturated | 
_ 5 | 300 | 0 | -175 | -119.7329 | 15.0601 | башабй | 
9 | 300 | 0 | -200 | -119.6301 | 14.9573 | башаейб —— | - - 
| 10 | O | SOk | -150 | -128.1796 | 23.5068 | Saturated __| | 
| 11 | 0 | 25k | -150 | -125.7621 | 21.0893 |  Saturaed_ | — | 
| 12 | 0 | 20k | -150 | -126.0112 | 21.3384 | Saturated | 
13 | 0 | 15К | - 1500 126167518922:00e A saturated ШШ 
| 14 | 0 | 10ок | -150 | -128.9777 | 24349) Saturated | 
| 15 | O | 8k | -150 | -129.6968 | 25.0040 |  Saturated | 
|. 16 | 0 Sk |-150 | -123.0963 | 18.425 | башабйб | | 
..17 | 0 | 8k | -75 |-118.3636 | 13.6908 |  Unsaturated | 
| 18 | о | 5% | -75 | 01677 | 54249 |. 
19 | 0) 2К OI A лы 
| 20 | 300 | 200” |-100 | -119.2771 | 146043 | AAA 
| 21 | 3004500 |-100|-1193445 | 1467217 | ОБ 
| 22 | 300 | 1000 | -100 | -104.6185 | 0.0543 |  Questionable | | 
_ 23 | 300 | 2К |-100 | - 96143 | 1494 | 7455 
24 | 300 | 4k | -100. -120.0087 | 15:36 | и 
| 25 | 500 | 6k | -100 | -120.3703 | 15769: EM ON 
_ 26 | 300 | 8k |-100|-1207338 | 160611 122 
27 | 300 Ок | 100 =121.0980 | 164259 AI 
728773009 12қ |2100: |2121:3796 | FOSA TN NE 

29 | 300 | ок | -100 | 12174791 "ИДР ae 
| 30 | 300 | 20k -100 | -1222401 |_ 175674) ш шиш 
231 | 3008| 7256721007] 212249732 TELAS н |ү 
| 32 | 3007] 50ке E1007) 21222592770 Ото | EA ИИ 
133 | 300 | 35 | -100 ЕРЛИУ БЛ 7 --. 
34 | 300 | 40к | -100 | -122.7332 | 18.0605 | Best Case This Series | | 


Table 25. ANSYS Series 2 Results 
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Series 3: 16.75 Hz 


oyo po | 0 -40622| - | Uncontrolled Case | l 
[1 | 30 | o | -20 | -136.3002 | -43220 |  AmplifiesSignal — | —— 
2 | 300 | 0 | -40 | -134.5450 | -6.0772_|  AmplifiesSignal | | 
__3 | 300 | 0 | -50 | -134.0553 | -6.5670_|  Amplifies Signal | У | 
[4 | 300 | 0 | 25 | -160.2686 | 19.6464 | Starting to Destabilize | У | 

5 | 300 | 0 | 50 |-142.6854 | 2.0631 | Negligible Effect_ | У | 
6 | 300 | 0 | 75 | -140.5569 | -0.0654 | Unstable to(-) limit | 
| 7 | 300 | 0 | 100 | -140.9973 | 03750 | Unstable | | 
[8 | 300 | 0 1125 | -141.0681 | 0.4459 |  Untbe | 
_9 | 300 | 0 | 150 | -141.1044 | 041) (ае | | 
_10 | 300 |_ о | 175 | -140.8921 | 02098) вые | | 
[11 | 300 | 0 | 200 | -141.1225 | 0.5003 | Unstable | 
{12 | 0 | 20k | 150 | -140.7954 | 0.1731 | Unstade | 
13 | 0 | 15 | 150 | -140.5506 | 0.0717 l 
| 14 | 0 | 10 | 150 | -140.5827 | 0.0395 | ае | 
| 15 | 0 | 5k | 150 |-141.0630 | 04408 | Unstable | 
_ 16 | 0 | 5k | 75 |-1414760 ^. 0.8537. | Unstable | 
[17 | 0 | 2k | 75 |-1414851 | 0.8628 | Starting to Destabilize | | 
| 18 | о | \ | 75 |-141.1150 | 0.4928 | Negligible Effect | _ 
| 19 | 0 | 500 | 75 | -140.8755 | 0.2533 | Negligible Effect | | 


Table 26. ANSYS Series 3 Results 
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Tnal 3: IGain=300, [[Gain=0, SG= -50. 


PCB Force Sensor and Controller Output 
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Figure 65. Series 3 Trial 3 Results and PSD 
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Trial 4: IGain-300, IIGain-0, SG= 25. 
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Figure 66. Series 3 Trial 4 Results and PSD 
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Trial 5: IGain-300, IIGain-0, SG- 50. 


PCB Force Sensor and Controller Output 
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Figure 67. Series 3 Trial 5 Results and PSD 
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Series 4: 11.30 and 11.75 Hz 
IGain-300, IIGain-0, SG= -75 


PO E A Ime A 
IAE EE ETICO AA EA 
о ОНГ к Me 
БЛ Б ee ЕОР Б г” 
E БО Es O PS E | и 
it Te we asi ar 
Rina a AA 721 
41 | 180 | -156.9814 | 14.8405 | 2024 
Boe РОО ОНИ 
SPA AA A O NI 
PAE PEA Жы ЖЫ” 
POMBO | E gi 


Table 27. ANSYS Series 4 Results 












185 


Amplitude (V) 


Trial 4: IGain=300, IIGain=0, SG= -50 Phase =120 degrees. 
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Figure 68. Series 4 Trial 31 Results and PSD 
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Power Spectrum Magnitude (dB) 


Series 5: Random Noise (0-200Hz) 
IGain=300, IIGain=0, SG= -75 
Peak Uncontrolled Amplitude = -199.2025dB 


Peak Controlled Amplitude = -205.2396dB 
Reduction = 6.0371dB 
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Figure 69. Series 5: Random Power Spectral Density 
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